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ABSTRACT 


Two  new  variational  principles,  of  particular 
interest  in  obtaining  information  about  excited  states, 
are  developed  and  applied  to  simple  systems. 

The  first  is  a  variational  principle  for  obtaining 
upper  and  lower  bounds  to  eigenvalues  of  the  time 
independent,  non-relativistic  SchrOdinger  equation. 

It  is  based  on  the  Temple-Kato  formulae  and  is  applied 
to  the  H*  ground  state  and  the  ground  and  excited  states 
of  the  He  atom.  A  comparison  is  made  between  this 
method  and  a  method  recently  proposed  by  Lowdin.  The 
methods  give  very  similar  results  but  the  present  method 
has  the  advantage  that  it  can  be  readily  extended  to 
more  complex  systems . 

The  second  is  a  variational  principle  which  results 
in  an  approximate  wave  function  for  the  ground  or  an 
excited  state  which  has  the  maximum  overlap  with  the 
true  eigenfunction.  The  method  involves  the  calcula¬ 
tion  of  E  =  <$|H[$>  and  A  =  <$|(H  -  E) 2 | $>  ,  but  is  free 
of  the  constraint  that  the  trial  function  $  must  remain 
orthogonal  to  all  states  of  the  same  symmetry  which  lie 
beneath  it.  One  must  know,  however,  an  approximation  to 
the  true  eigenvalue,  but  a  discussion  is  given  on  how 
one  might  gain  the  latter  information,  lacking  knowledge 
of  the  spectrum,  by  repeated  application  of  the  method. 


■ 


* 


The  method  is  applied  to  some  excited  states  of  the 
He  atom.  It  is  shown,  on  the  basis  of  physical  properties 
calculated  from  the  resultant  wave  functions,  that  this 
method  is  superior  to  the  familiar  variational  principle 
used  in  most  quantum  chemical  calculations.  As  a  result 
of  this,  it  is  hoped  that  this  new  method  may  prove 
itself  to  be  of  particular  value  in  the  future  in  treating 
problems  of  chemical  interest  and  determining  accurate 
and  tractable  wave  functions. 
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I.  General  Theory  of  Upper  and  Lower  Bounds  to  Eigenvalues 

1.1  Introduction 


In  the  application  of  the  Quantum  Theory  to  systems 
of  chemical  interest,  one  is  usually  attempting  to  find 
an  approximation,  $,  to  the  true  wave  function,  y,  of  the 
time  independent  non- relativistic  Schrodinger  equation. 
That  is,  one  is  trying  to  solve  the  equation 

Hy  =  Wy  (1) 


and  determine  the  energy,  W,  of  a  system  and  the  wave 
function  associated  with  W;  but  one  must  in  general 
compromise  with  an  approximate  $  and  an  approximation 

to  W. 


For  the  ground  state,  the  variational  principle  states 


that 


E  =  < $  I H I $>  >  W. 


(2) 


This  gives  an  upper  bound  to  the  true  energy  W,  but  it 
gives  no  estimation  of  the  error  when  the  true  W  is  un¬ 
known,  since  there  is  then  no  lower  bound  to  the  energy. 
It  is  of  the  utmost  importance  for  the  prediction  of 
physical  phenomena  that  there  exist  means  to  obtain  both 
upper  and  lower  bounds  to  eigenvalues,  not  only  of  the 
ground  state,  but  also  excited  states.  It  is  the  purpose 
of  this  work  to  1)  investigate  various  criteria  for 
obtaining  upper  and  lower  bounds  and  to  formulate  a 


■ 


method  based  on  a  variational  principle  for  the  determina¬ 
tion  of  these  bounds,  and  2)  develop  a  new  variational 
principle  for  excited  states. 

In  section  1.2  the  formulae  which  have  been  proposed 
to  determine  an  upper  or  lower  bound  from  a  knowledge  of 
a  variat ionally  determined  approximate  wave  function  4^. 
will  be  reviewed. 

This  will  lead  in  section  1.3  to  the  question:  what 
is  the  $  which  serves  as  the  best  approximation  to  ¥? 

After  reviewing  criteria  for  the  quality  of  a  trial  function, 
a  method  will  be  given  to  obtain  the  best  i.e.,  the 

best  approximation  to 

After  having  obtained  a  method  for  determining  the 
best  in  section  1.4, the  question  of  the  best  bounds 

is  taken  up.  It  is  discovered  that  the  combination  of  the 
best  ^  determined  variationally  by  the  approach  given 
in  section  1.3  when  combined  with  the  bound  formulae 
discussed  in  section  1.2  does  not  give  the  best  upper 
and  lower  bounds.  A  method  is  then  presented  to  obtain 
the  best  bounds.  It  is  noted  there  that  the  formulae 
for  the  best  bounds  and  the  formula  for  the  best 
discussed  in  section  1.3,  have  the  same  form. 

In  section  1.5  consideration  of  the  work  in  the 


literature  is  given.  It  is  shown  that  the  methods  of 
sections  1.3  and  1.4  provide  a  good  basis  for  an  analysis 
of  the  work  in  the  literature.  Also  the  literature 
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supports  the  contentions  made  in  sections  1.3  and  1.4. 

In  section  1.6  alternative  methods  to  the  calcula¬ 
tion  of  lower  bounds  are  considered.  All  of  the  work 
considered  in  sections  1.2  -  1 . 5,  and  indeed  all  the  work 
which  is  presented  in  Chapters  II  and  III,  are  based  on 
a  variational  principle.  In  section  I. 6, however,  methods 
based  on  a  perturbation  approach  are  outlined  for 
completeness . 


I . 2  Formulae  for  Bounds 

The  first  mention  in  the  literature  of  a  method  for 
determining  a  lower  bound  to  eigenvalues  is  in  the  work 
of  G.  Temple  (1).  He  proposed  that,  for  the  ground  state, 
one  should  obtain  a  lower  bound  by  applying  the  formula 

B  =  E  +  tr- - —  (3) 

h  -  y  v  ' 

where  B  is  the  lower  bound  and 

E  =  <$|H|$>,  A  =  <$|(H  -  E) 2 | $>  (4) 

and  where  y  is  an  approximation  to  the  eigenvalue  of  the 
first  excited  state.  $  is  an  approximate  eigenfunction 
which  satisfies  Eq.  (2).  Thus  Eq.  (2)  gives  an  upper 
bound  and  Eq.  C3)  a  lower  bound  to  the  exact  eigenvalue. 
Later  D.  H.  Weinstein  C2)  proposed  another  method 
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for  lower  bounds  which  could  also  be  used  for  excited 
states.  According  to  Weinstein  there  is  some  energy 
level  Wj,  which  satisfies  the  relation 

E  +  /a"  >_  W^.  >  E  -  /a"  (5) 

where  E  and  a  are  as  defined  in  Eq .  (4).  A  may  also  be 
written  as 

A  =  D  -  E/_ ,  where  D  =  <H$|H$>  (6) 


If  $  is  expanded  in  the  complete  set  of  eigenfunctions 


¥  - ,  then 


$  =  [  a  'F 
u  n  n 
n 


(7) 


and  thus 


E  =  \  a* a  W  ,  D  =  \  a*a  W2  and  \  a*a  =  1.  (8) 

‘■‘nnn’  ^  n  n  n  Lnn  K  J 

n  n  n 


In  what  follows  it  will  be  assumed  that  $  is  a  real 

r\ 

function  and  hence  a*a  will  be  written  as  a  .  From 

n  n  n 

Eqs .  (6)  and  (8)  one  obtains 


A  =  D 


I  a2W2 
L  n  n 
n 


2E  l  a2W  +  E2  l  a2  =  I  a2  (W  -  E) 

Lnn  ^  n  Ln^n  J 

n  n  n 

(9) 


The  argument  then  is  that  there  is  some  energy  level  W^ 
which  is  at  least  as  near  as  any  other,  i.e.,  for  which 


(Wk  -  E)2  1  CWn  -  E) 2 


CIO) 


A.  >  (W,  -  E)2  l  a2  or  A  >  (Wfc  -  E)  2.  (11) 

n 


Hence 


. 
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There  are  now  two  cases  to  consider: 

W^  >_  E  and  <  E. 

In  the  first  case 

/A  >  Wk  -  E  and  thus  E  +  /a"  >_  >_  E  (12) 

and  in  the  second 

/a  >_  E  -  and  E  >  >_  E  -  /a.  (13) 

Combining  Eqs .  (12)  and  (13),  Weinstein’s  statement, 

Eq .  (5),  is  obtained.  This  method  however  has  some 
undesirable  characteristics  which  will  be  pointed  out 
later . 

Several  other  methods  for  calculating  lower  bounds 
based  on  the  use  of  A  and  E  were  proposed,  but  the  work 
of  Kato  (3),  placed  all  of  them  on  a  common  basis.  Indeed, 
he  showed  that  the  Temple  formula  was  the  best  lower 
bound  which  could  be  obtained  with  a  knowledge  of  A, 

E  and  an  approximation  to  the  energy  of  the  first  excited 
state.  He  moreover  extended  the  method  to  excited  states, 
to  yield  upper  bounds  as  well  as  lower  bounds.  Kato’s 
theorem  will  now  be  stated. 

Kato's  Theorem:  If  it  is  assumed  that  the  interval 
(E^  ,  Ej.+-^)  contains  no  point  of  the  spectrum  other  than 

a  non- degenerate  eigenvalue,  and  that,  for  a  trial  func¬ 
tion  $,  E  and  A  are  calculated  and  satisfy  the  inequality 
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A  <  (E  -  E.^KE.^  -  E)  C14) 

th  en  the  interval  certainly  contains  one  and  only  one 
eigenvalue  VT  such  that 


E  + 


E  -  Ei+i 


=  B*  <  W.  <  BU  =  E  + 
—  1  — 


E  -  E.  , 
l-l 


(15) 


In  general  Eq.  (15)  cannot  be  replaced  by  a  more 
precise  formula  if  only  the  quantities  E,  A,  E^  ^  and  E^  +  ^ 
are  known,  whatever  these  quantities  may  be.  The  proof 
is  as  follows : 

One  may  write  the  Hamiltonian  in  terms  of  its  eigen¬ 
values  and  eigenfunctions  as 


& 

h  =  y  w.  y. 

v  i  i  i 

where  it  is  assumed  here,  for  simplicity,  that  the 
is  entirely  discrete. 

It  is  assumed  that  for  a  normalized  vector  $ 
domain  of  H  the  following  quantities  are  calculate 


(16) 

spectrum 


in  the 
d 


<$|H|$>  =  E,  <$|(H  -  E) 2 | $>  =  A  .  (17) 

Lemma  1 :  For  every  E^_^  such  that  E^_^  <  E  (E^^  -°°  is 
permitted),  the  interval  (E^_1,  Bu]  contains  a  point  of 
the  spectrum,  where  Bu  is  defined  by 


(18) 


- 
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Proof  of  Lemma  1:  If  the  interval  CE^  ^ for 

y  >  Ej:_2>  contains  no  point  of  the  spectrum,  the  quadratic 

expression 

(E  -  Ei_1)(E  -  y)  =  E2  -  (Ei_1  +  y)E  +  E^  ^  y 

(19) 

cannot  take  on  a  negative  value  as  E  varies  over  the 
spectrum.  The  same  will  be  true  if  E  is  replaced  by 
and  then  Eq.  (19)  is  multiplied  by  < $  |  | $> ,  and 

this  expression  is  summed  over  i.  Since 


and 


it  follows  that 


or 


<<F  .  1  $>  = 
1  1 

<$ | H2 | $>  =  A  +  E2 

(20) 

<$ 1 >w. 

1  l  l 

<y . 1 $>  =  E 

l 1 

(21) 

-  Cq.! 

+  y )  E  +  E^_i  y  >_  0 

(22) 

<  E  + 

A  =  Bu. 

(23) 

E  -  E. 

l-l 


This  means  that  the  interval  (E^_^,  y)  contains  a  point 
of  the  spectrum  if  y  >  Bu,  and  since  the  spectrum  is  a 
closed  set,  the  same  is  true  of  the  interval  (E^_^,  Bu]  . 


Lemma  2 :  For  every  E^  +  ^  such  that  E^+-^  >  E,  the 

a 

interval  [B  ,  E^+^)  contains  a  point  of  the  spectrum, 
where  B  is  defined  by 


=  E  +  - — 


E  -  E , 
1  +  1 


C24) 


V 


• 


' 
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This  proof  follows  the  one  given  for  Lemma  1.  Now, 

combining  the  two  lemmas, it  is  clear  that  each  of  the 
**  u  £ 

intervals  (E^_^,  B  ]  and  [B  ,  E^  +  ^)  contains  a  point  of 
the  spectrum.  Since  equation  (14)  really  means  that 
BU  <  E.  ,  and  B5,  >  E,  ,  ,  it  must  be  concluded  that  [b\  Bu] 
contains  a  point  of  the  spectrum,  . 

In  the  above  theorem  B  and  B  are  the  lower  and  upper 

bounds  respectively,  and  the  E  are  approximations  to  the 

th 

eigenvalues  of  the  yL  state. 

Hereafter  we  will  consider  in  detail  only  Weinstein's 
method  and  Temple’s  method  as  extended  by  Kato  as  a  means 
of  calculating  bounds,  since  the  former  gives  the  best 
bound  obtainable  with  the  information  E  and  A,  and  the 
latter  gives  the  best  bound  obtainable  from  E,  A  and  E^ . 
There  are  other  methods,  which  will  be  mentioned  briefly 
later,  which  give  lower  bounds  without  the  knowledge  of 
A.  These  methods  avoid  the  calculation  of  A,  which  involves 
quite  difficult  integrals,  but  introduce  complications 
of  their  .own. 

It  should  be  noted  that  the  formulae  which  will  be 
considered,  namely  Eqs .  (5)  and  (15),  assume  that  one 
knows  an  approximate  wave  function  to  the  k^*1  state, 
in  order  to  obtain  a  bound  to  the  eigenvalue  of 
that  state.  This  obviously  presents  the  problem  of 
determining  such  an  approximate  wave  function  for  the  k^^1 
state.  This  problem  will  now  be  considered  in  detail. 

\ 


. 


c  .«•  t  n  ;,  n*  *  (  ,  J  htn  [  fi  ,  ,  3)  t 


’ 


‘ 

1  « 
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1.3  Determination  of  . 

k 


For  the  exact  wave  function  y  of  an  eigenstate  the 
following  equations  hold, 


A 

3? 

•-€ 

V 

II 

(25) 

<T |H2  -  E2 | =  0 

(26) 

and  if  y  is  an  eigenstate  with  E  f  E',  where  E'  is  the 
energy  of  another  state  of  the  same  symmetry  with  the 
corresponding  eigenfunction  y T ,  then 

<Y | >  =  0  .  (27) 

For  a  trial  function  $,  one  has  E  =  <$|H|$>  >  W  ,  where 
WQ  is  the  true  ground  state  energy.  Indeed  if  $  were 
capable  of  attaining  the  exact  form  of  y  ,  the  exact  ground 
state  eigenfunction,  one  would  have  E  =  Wq  =  <^0 I H | y  > . 

Likewise  for  a  calculation  of  an  excited  state,  if 
the  trial  function  $  is  capable  of  attaining  the  form 
of  under  a  variation,  then  a  stationary  state  is 
obtained  with  E  =  I H  |  .  In  both  cases  it  is 

to  be  noted  that  the  true  eigenvalues  are  obtained  by 
finding  the  stationary  value  from  6E  =  0  for  the  case  of 
a  trial  function  which  can  become  the  eigenfunction  under 
variation . 

For  an  excited  state  of  the  same  symmetry  as  the 
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ground  state  ,  5E  =  s<$|H|$>  =  0  is  a  necessary  condition 
for  finding  a  stationary  value  E  which  approximates  W^; 
but  it  is  not  sufficient.  One  must  also  take  account 
of  Eq.  C 2 7 )  .  Thus  with  Eqs .  C 2 5 )  and  (27)  using  $  in 
place  of  y  one  has  necessary  and  sufficient  conditions 
so  that  $  +  5<$|H|$>  =  0  and  <$|t  >  =  0.  But  the 

condition  <$|$q>  =  0  is  rather  restrictive,  especially 
if  is  not  a  very  good  approximation  to  y  .  Hence 
it  is  desirable  to  find  another  approach  for  excited  states 

Now  Weinstein  (2)  has  proposed  that  for  an  excited 
state  for  which  one  has  a  trial  function  $,  which  cannot 
obtain  the  form  of  an  exact  wave  function,  one  can  obtain 
an  approximation  to  some  4^  by  minimizing  A,  where  A  is 
given  by  Eq .  (4) . 

A  simple  theorem  states  that:  if  there  exists  a  $ 
such  that  H$  =  W^.$  (hence  $  =  y^.)  ,  then  A  =  0.  As  a 
corollary  it  may  be  stated  that  as  $-*■  A  0 .  Hence 
$  -»■  implies  that  A  -*■  0 .  It  is  on  this  basis  that 
Weinstein  suggested  that  as  a  ^  min.,  $  -*  It  will 

be  shown  below  that  if  one  measures  the  approach  of  $  to 
'i'k  by  the  increase  of  the  value  of  the  integral  a^  = 
the  minimization  of  A  does  not  at  all  give  the  best 
approximation  to  a  y-^  that  can  be  obtained  from  a  given 
trial  function. 

There  are  criteria  for  judging  the  goodness  of 
$  as  an  approximation  to  other  than  that  of  the 
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maximum  a^,  but  a^  provides  for  the  best  overall  con¬ 
vergence  of  $  to  ^  •  Obtaining  the  maximum  a^  is 
equivalent  to  obtaining  the  smallest  root-mean-square 
deviation  of  $  from  .  As  will  be  seen  below,  the 
other  criteria  are  useful  for  obtaining  good  approximations 
to  4^  in  certain  regions  of  configuration  space,  but 
since  aj,  provides  a  criterion  for  judging  the  goodness  of 
$  as  an  approximation  to  4^  over  the  entire  configuration 
space,  it  seems  that  this  criterion  should  be  the  one 
used  to  judge  the  goodness  of  $  . 

Criteria  of  Goodness  for  Approximate  Wave  Functions 


In  this  section,  a  generalization  is  presented  of 
an  analysis  first  given  by  H.  M.  James  and  A.  S.  Coolidge 
(4)  for  the  ground  state  of  an  electronic  system.  Here 
the  method  will  be  extended  to  apply  to  excited  states  as 

well . 

Let  $  be  an  approximate  normalized  wave  function  and 
E  =  <$|H|$>  its  associated  energy.  Also  let 


Now 


$  =  T  a. 4'.  =  7  a.  4f.  +  a^v,  . 

i  1  1  ifk  1  1  k  k 


O  -  a^iJ  = 


vtW  -  >  a.  ¥• 

K  K  i£k  1  1 


(28) 

(29) 


will  give  the  deviation  of  $  from  the  exact  function  4^ 
A  deviation  function  can  then  be  defined, 


- 


. 
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*x  -  u  -  ak)  2c*  -  akV 


(30) 


where  a,  =  Hence  one  may  write 


4  =  Vk  +  axtx 


(31) 


2 

where  a  =  (1  -  a J)  measures  the  amount  of  the  deviation 

X  K 

function  4  which  appears  in  $.  As  criteria  of  the  in- 
accuracy  of  $  there  are:  Q,  the  root -mean- square  error 


m  $ 


Q  =  <  $ 


Y,  $ 

k ' 


v 


(32) 


e,  the  energy  error 


e  =  E  -  W, 


(33) 


and  A  ,  the  root -mean-square  local  energy  deviation 


1  1 

A2  =  <$! (H  -  E)  2  I  $>  2  . 


(34) 


One  may  also  define  the  quantities 


and 


(35) 


If  the  inequality 

<$x| [ CH  -  E)  -  CEX  -  E)]2|$x>  >  0 

is  considered,  and  the  integral  is  expanded  in  terms  of 
the  above  defined  quantities,  one  finds 


■ 


-13- 


<$  |H2|0  >  -  2E  E  +  E2  -  E2  t  2E  E  -  E2  >  0, 

X  X  X  X  X  — 

Thus  , 

<$  |H2U  >  -  2E  E  +  E2  >  E2  -  2E  E  +  E2. 
x1  1  x  x  -  x  x 

But  the  left  hand  side  is  merely  A  ,  hence 

7  x’ 

Ax  i  (Ex  -  E)2.  (36) 

From  the  definitions  of  e  and  e  given  above,  it  can  be 
shown  that  (E  -  E) 2  =  (e  -  e)2,  hence  A  >_  (e  -  e)2. 

X  X  XX 

Thus  one  may  define  a  quantity  K2 , 


K2  = 


(e  -  e) 
v  X  J 


>  1  . 


(37) 


Now  substituting  Eq.  (31)  into  Eqs .  (32)  -  (34), it 
follows  that  1 

Q2  =  2 [1  -  ak]  =  2 [1  -  (1  -  a2)2]  (38) 


£ 


a2  e 

X  X 


(39) 


(40) 


From  the  first  equality  of  Eq .  (38)  one  finds  the 
justification  of  the  statement  made  in  the  last  section, 
that  obtaining  the  maximum  a^  is  equivalent  to  obtaining 
the  smallest  root-mean-square  deviation. 

Eliminating  a  between  Eqs.  (38)  and  (39)  one  finds 

X 


(41) 


it  ■  j  K  ii  1 3m 


■ 
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or  if  $  is  a  fairly  good  approximation 
£  is  small  when  compared  to  | ^ 
then 


to  ^  (i . e . ,  if 
and  |Wk  -  Wk+1|), 


(42) 


Eliminating  a  between  Eqs.  (39)  and  (40)  one  obtains 

X 


A  =  s  [e  +  (Ax  -  e 2  3 /e  x] 


(43) 


and  assuming  again  that  $  is  a  fairly  good  approximation 
to  it  follows  that 


A/A  -  e/e 
x  x 


or  .  A  -  ee  Kz . 

x 


(44) 


Now,  making  use  of  Eqs.  (42)  and  (44),  one  obtains  the 
result 

o  £  ^  0 

(45) 


Q2  =  (— )K2. 

A 


This  equation  interrelates  the  various  criteria  considered. 
Again  it  should  be  emphasized  that  theoretically  the  criterion 
which  should  be  used  to  judge  the  overall  quality  of  a 
wave  function  is  the  criterion  of  minimum  Q2  (i.e., 
maximum  a^.)  . 

Before  proceeding  to  investigate  Eq.  (45),  each  of 
the  quantities  in  this  equation  must  be  studied  in  order 
to  understand  its  consequences  and  enable  one  to  proceed 
to  a  method  for  determining  the  best  $  to  approximate 

Assume  first  that  one  has  a  function  $.  One  may 
then,  according  to  Eqs.  (32)  -  (34),  judge  its  accuracy. 
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It  should  be  noted  that  in  applying  these  formulae  directly 
that  Eq .  (32)  requires  a  knowledge  of  the  true  wave  func¬ 
tion.  This  is  not  available,  of  course,  except  for  the 
very  simplest  systems;  thus  Eq .  (32)  is  impractical  from 
this  point  of  view. 

In  order  to  use  Eq .  (33)  one  must  know  the  true 
eigenvalue  ;  however,  for  this  quantity  the  experimentally 
determined  energy  for  this  state  may  be  used,  if  it  has 
been  measured.  Thus  this  method  of  judging  accuracy 
usually  will  be  available. 

Finally,  it  should  be  noted  that  Eq .  (34)  makes  no 
use  of  information  other  than  that  contained  in  $. 

Although  this  seems  to  be  a  favorable  aspect  of  Eq .  (34) 
it  must  be  realized  that  one  is  only  finding  the  local 
energy  deviation  on  the  basis  of  $  without  regard  to  any 
external  standard. 

The  three  criteria  have  now  been  considered  assuming 
that  a  $  had  somehow  been  obtained.  A  more  important 
question  is  whether  the  criteria  can  be  used  to  determine 
a  $.  This  is  obviously  the  only  way  one  can  obtain  a 
$  which  approximates  a  because  if  one  does  not  have 

criteria  to  guide  one  in  obtaining  a  $,  how  then  does 
one  obtain  such  a  function? 

From  this  point  of  view  Eq.  (32)  is  useless;  because 
if  is  known,  why  should  one  want  to  obtain  an  approxima¬ 
tion  to  it?  Eq.  (33)  can  be  used  to  obtain  a  4,  if  one 


- 
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assumes  that  $  is  the  function  which  minimizes  e.  That 
is ,  a  $  must  be  found  such  that  Se  =  0.  For  the  ground 
state  this  is  nothing  more  than  the  familiar  variation¬ 
al  principle.  Eq.  (33)  is,  however,  incapable  of  giving 
a  $  which  is  an  approximation  to  an  excited  state  which 
has  the  same  symmetry  as  the  ground  state. 

Considering  Eq.  (34)  one  must  find  the  $  which 
minimizes  A.  it  is  the  only  one  of  the  criteria  which, 
taken  alone,  can  yield  a  $  for  the  ground  state  or  any 
excited  state.  It  is  for  this  reason  that  it  has  been 
suggested  by  Weinstein  as  a  method  for  determining 
approximations  to  .  Referring  again  to  Eq .  (45)  however, 

it  is  obvious  that  Weinstein's  idea  of  minimizing  A  does 
not  give  the  minimum  Q2  and  thus  this  method  cannot  give 
the  best  overall  approximation  to  .  But  looking  at 
each  of  the  three  criteria  individually  and  demanding  that 
Q  be  a  minimum,  an  impasse  seems  to  be  reached.  It  will 
be  shown  shortly  that  this  impasse  is  resolved  by  Eq .  (45). 
But  first  brief  consideration  will  be  given  to  the  types 
of  error  in  the  wave  function  which  are  reduced  by  minimiza¬ 
tion  of  A  and  the  minimization  of  E;  here  the  discussion 
given  by  James  and  Coolidge  (4)  is  important. 

In  minimizing  A,  the  "local  energy  error"  is  minimized 
(4,  5).  That  is,  the  errors  in  the  wave  function  associated 
with  the  closeness  of  the  electrons  will  be  minimized. 

Thus  the  wave  function  which  results  from  this  type  of 
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minimization  would  be  expected  to  giye  good  values  for 
properties  which  describe  phenomena  where  the  particles 
are  close  together. 

In  minimizing  E,  which  can  be  done  only  for  the  ground 
state  or  lowest  state  of  a  given  symmetry,  "long  range 
errors”  are  minimized  (4,  5,  6).  Thus  a  wave  function 
determined  on  this  basis  should  be  good  for  describing 
properties  such  as  the  radial  density  or  diamagnetic 
susceptibility  where  long  range  errors  are  important. 

Thus  by  minimizing  E  or  A,  a  wave  function  is  obtained 
which  has  minimum  error  in  certain  regions  of  configuration 
space,  rather  than  the  minimum  error  over  the  whole  of 
configuration  space. 

Now  returning  to  Eq.  (45)  one  may  ask  how  one 
might  obtain  a  wave  function  for  an  excited  state  which 
satisfies  the  criteria  of  minimum  Q2  .  From  the  deriva¬ 
tion  of  Eq.  (45)  it  was  found  that  K2  1,  according 
to  Eq .  (37).  Thus,  if  the  right  hand  side  of  Eq .  (45) 
is  to  be  minimized,  the  smallest  value  K2  can  take  is  1. 
Thus,  in  order  to  obtain  the  minimum  Q2 ,  the  quantity 
e 2 / A  should  be  minimized.  This  method  allows  the 
determination  of  an  excited  state  wave  function  which 
minimizes  Q2  without  a  knowledge  of  the  exact  This 

seems  to  be  the  first  time  that  this  procedure  has  been 
proposed  as  a  method  for  obtaining  approximations  to 
excited  state  wave  functions. 


■ 
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On  the  basis  of  the  new  method  just  proposed  above, 
it  is  now  known  how  to  obtain  a  $  which  is  the  best 
approximation  to  a  It  should  be  noted  that  the 

obtained  by  this  method  for  the  case  of  a  ground  state 
need  not  give  as  good  an  approximation  to  W  as  could 
be  obtained  by  minimizing  E;  likewise,  it  need  not  give 
as  small  a  value  for  A  as  would  be  obtained  by  minimizing 
A.  However,  unlike  the  minimization  of  A  or  E,  it  does 
give  the  minimum  Q2  which  is  the  chosen  criterion  for 
the  best  . 

I . 4  The  Determination  of  the  Best  Bounds 

In  the  last  section  it  was  found  that  one  can  obtain 
the  best  by  minimizing  e2/A.  Hence  with  this  wave 
function  one  may  return  to  the  formulae  of  section  1.2 
for  upper  and  lower  bounds.  In  particular  one  is  interested 
in  Eqs .  (5)  and  (15).  If  is  substituted  into  Eq .  (15) 
along  with  the  E  ,  the  best  bounds  to  the  kL  state 
will  be  obtained  that  ^  can  give.  If  there  is  no 
knowledge  of  the  E  ,  then  one  must  settle  for  the  bounds 
from  Eq .  (5)  . 

Thus  the  best  bounds  available  for  the  function 
have  been  obtained.  But,  have  the  best  bounds  to  the 
k^*1  state  which  can  be  obtained  from  the  form  of  the 
trial  function  $  been  determined?  The  answer  is  no, 
not  necessarily.  In  the  last  section  it  was  found  that 
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the  best  function  judged  by  A  minimization  is  not 
the  same  as  the  best  function  as  obtained  by  e  minimiza¬ 
tion  of  e2/a;  similarly  the  bounds  based  on  obtained 

by  Eq .  (45),  may  not  be  as  good  as  the  ones  based  on 
another  method.  The  method  to  be  used  to  determine  these 
better  bounds  will  be  outlined  now. 

First,  consider  Eq .  (5);  for  the  lower  bound  one  has 

B£  =  E  -  /a  (46) 

and  for  the  upper  bound 

Bu  =  E  +  /a  .  (47) 

If  one  wishes  to  obtain  the  best  bounds  which  can  be 
given  by  a  function  of  the  form  $ ,  one  should  maximize 
Eq .  (46)  with  respect  to  $  and  obtain  a  $  *  and  then 

minimize  Eq.  (47)  with  respect  to  $  and  thus  obtain  a 
$  .  If  $  and  $  are  determined  with  regard  to 

the  k1'*1  state  they  may  be  labeled  as  $™ax*  and  ‘  . 

The  relationship  which  exists  between  $™ax*,  ’  and 

is  at  best  obscure.  Thus  it  seems  that  one  is  given 
the  choice  of  obtaining  the  best  anc*  calculating 
bounds  from  it  which  will  not  necessarily  be  the  best 
bounds,  or  calculating  the  best  bounds  and  obtaining 
thereby  $^ax*  and  which  may  not  necessarily  be 

good  approximations  to 

At  this  point  the  choice  of  trying  to  determine 
the  best  bounds  to  a  given  state  will  be  made,  rather 


. 
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than  trying  to  determine  the  best  What  is  clear 

from  the  preceeding  discussion  is  that  the  method 

proposed  by  Weinstein,  i.e.,  the  minimization  of  A, 

yields  neither  the  best  nor  the  best  bounds.  This 

point  will  be  returned  to  again  later. 

Having  chosen  this  course,  it  is  necessary  now  to 

consider  the  Temple-Kato  formula,  Eq .  (15).  Here,  as 

with  the  Weinstein  formula  considered  above,  one  wishes 

to  maximize  the  lower  bound  and  minimize  the  upper 

bound.  Hence  one  again  will  obtain  a  $^ax*  and  a 

With  Eq .  (15),  if  a  maximum  and  minimum  can  be 

found  at  all,  it  will  be  guaranteed  by  Kato's  theorem 

that  the  corresponding  bounds  B£  and  Bu  will  be  bounds 
*th 

to  the  kL  state  which  is  under  consideration.  This 
is  in  contrast  to  the  method  based  on  the  maximization 
of  Eq.  (46)  and  the  minimization  of  Eq.  (47).  For  although 

£  U 

one  will  obtain  in  this  manner  the  best  B  and  B 
obtainable  from  Weinstein's  formula,  one  is  not  insured 
that  they  will  be  the  best  bounds  to  the  state  under 
consideration;  they  may  in  fact  be  the  best  bounds  to 
an  entirely  different  state.  In  order  to  get  the  best 
bounds  to  the  desired  state,  one  must  make  a  very 
judicious  choice  of  the  trial  function  $ .  This  is  one 
inconvenience  of  a  theory  based  on  the  Weinstein  formula. 

It  is  to  be  noted  that  the  bounds  obtained  from 
maximizing  B2,  and  minimizing  Bu  as  given  by  Eq.  (15)  are 
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in  any  event  always  better  than  those  obtained  from 
maximizing  Eq.  (46)  and  minimizing  Eq.  (47),  (3).  Of 
course  in  operating  with  Eq.  CIS)  one  must  always  have 
a  knowledge  of  the  E  . 

It  is  now  proposed  that  the  method  outlined  above, 
which  is  based  on  Eq .  (15),  will  determine  the  best 

bounds  which  can  be  obtained  from  a  knowledge  of  E  and 
A.  Now  proceeding  from  the  outline  given  above,  this 
new  proposal  will  be  considered  in  detail. 

Rewriting  Eq.  (15)  in  the  following  form 

A 

B  =  E  +  - —  (48) 

E  -  E 

Y 

where  B  is  an  upper  or  lower  bound  depending  upon  whether 
Yisk-lork+1,  and  demanding  the  variation  in  B 
should  vanish  for  the  bounds,  the  following  is  obtained 

<5B  =  0  =  SE  +  — — - - .  (49) 

E  -  E  (E  -  E  )2 
Y  Y 

This  may  be  rewritten  as 

SB  =  SE  +  [ - ]  6 A  =  0  (50) 

(E  -  E  )2  -  A 
Y 

which  states  that  the  trial  function  $  which  makes  the 
variation  in  B  vanish  will  serve  as  the  function  to 
obtain  the  best  bound  when  it  is  used  in  Eq.  (48) . 

In  particular,  it  is  assumed  that  one  is  interested 
in  obtaining  upper  and  lower  bounds  to  the  eigenvalue  of 


■ 
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the  k**1  state.  For  the  upper  bound  of  the  k1"^  state 

one  would  need  Ek_^,  an  approximation  to  the  eigenvalue 

til  ~ 

of  the  (k  -  1)  state.  With  Ek_^  and  a  4,  one  would 
use  Eq.  (50)  to  determine  the  $^in * ,  which  gives  the  best 
upper  bound  to  W^ .  The  $  which  satisfies  Eq .  (50)  in 
this  case,  $kin‘,  must  be  such  as  to  yield  the  minimum 
value  of  Eq .  (48) . 

For  the  lower  bound  one  would  use  E,  , ,  and  determine 

k+1 

by  use  of  Eq.  (50)  the  $^ax*  which  would  yield  the  best 

lower  bound.  This  function  would  yield  the  best  lower 

bound  obtainable  from  a  function  of  type  $ ,  when  it  was 

substituted  into  Eq.  (48) . 

Thus,  so  long  as  one  knows  an  Ek+-^  and  an  E^  ^ , 

one  is  insured  by  Kato’s  theorem  that  when  Eq.  (50) 

t  ll 

is  used  one  will  determine  bounds  to  the  kL  state; 

Eq .  (50)  also  guarantees  that  the  bounds  are  the  best 

obtainable  from  a  function  of  the  form  $.  It  should 
be  noted  that  from  Kato's  theorem  the  inequality 
A  <  (E  -  E^  ^)  (Ej,  +  ^  "  E)  must  be  satisfied.  This  does 
not  however  imply  that  A  must  bracket  W^,  i.e., 
E-/a<W-,  <E  +  /A,  as  seems  to  have  been  assumed  in 
the  literature  in  the  past  (7,  8,  9).  On  the  contrary 

m 

A  as  determined  by  a  ^  *  may  be  such  that 

E  -  /a  <_  Wk+1  1  E  +  /a,  and  yet  this  A  and  $>kax’  will 

th 

serve  to  give  a  lower  bound  to  the  k  state.  The  only 
restriction  on  A  is  that  it  must  satisfy  the  inequality 
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Eq.  (14) ,  in  Kato's  theorem. 

Now  Eqs .  (46)  and  C4 7 )  may  be  written  as 

B  =  E  ±  /A  (51) 

where  B  is  an  upper  or  lower  bound  depending  upon  whether 
the  plus  or  minus  sign  is  chosen.  There  is  an  interesting 
relationship  among  Eqs.  (45),  (48)  and  (51)  which  is 
worthwhile  developing. 

Considering  Eq.  (51)  first,  if  the  variation  in 
B  is  set  to  zero  one  obtains 

5B  =  6E  +  [ — — ]  6 A  =  0  (52) 

2/A 

This  equation  is  to  be  interpreted  in  the  same  manner  as 
Eq .  (50);  note  also  the  similarity  of  the  two  equations. 
They  differ  only  in  the  coefficient  of  6A. 

Now,  considering  Eq .  (45),  the  variation  in  Q2 

should  vanish.  Hence  one  writes,  considering  the 
arguments  put  forward  before  on  the  minimization  of  Q2 , 

6Q2  =  <5E  -  [ — - ]6A  =  0  (53) 

which  has  the  same  form  as  Eqs.  (52)  and  (50).  The  only 
difference  is  again  the  coefficient  of  <5A.  Thus  one 
may  write  all  three  of  these  equations  in  the  common 
form 

6E  +  Co  6  A  =  0  (54) 

where  the  coefficient  to  depends  upon  whether  one  wants 
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a)  the  best  bounds  obtainable  from  a  $  with  knowledge 
of  E  and  A  and  E  ,  Eq .  (50) 

b)  the  best  bounds  available  from  a  $  with  a  knowledge 
of  only  E  and  A,  Eq .  (52) 

c)  the  best  wave  function  obtainable  with  the  knowledge 

of  only  Eq.  (53). 

Thus  all  of  the  methods  which  have  been  developed 
can  be  summarized  by  Eq .  (54).  Whether  there  is  any 
deeper  significance  to  the  fact  that  these  three  methods 
can  be  represented  by  a  common  form  is  at  the  moment 
not  clear.  But  it  should  certainly  be  interesting  to 
pursue  this  line  of  thought  further. 

Before  leaving  this  section  on  the  determination 
of  the  best  bounds  it  may  be  of  profit  to  mention  a  draw¬ 
back  in  the  use  of  Eqs .  (48)  and  (50),  since  these  are 
the  equations  which  will  be  used  to  develop  a  working 
scheme,  based  on  orbital  theory,  in  the  next  chapter. 

The  drawback,  of  course,  is  that  one  has  to  know  the 
Ey .  The  better  they  are  known,  the  better  the  bounds 
will  be.  The  best  values  are,  of  course, E^+1  =  W^+-^ 
and  E^  ^  =  ^k-lJ  anc*  t^lus »  if  there  is  experimental 

information  available,  it  could  be  used  for  the  E  . 

Y 

However,  it  may  be  that  no  experimental  information  is 

available  in  a  particular  case;  how  then  could  one 

obtain  the  E  ? 

Y 

The  E  could  perhaps  be  determined  a  priori  in  the 


. 
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following  way.  Suppose,  for  a  definite  example,  that 
one  wishes  to  determine  the  upper  and  lower  bounds  of 
the  second  excited  state  of  a  given  symmetry  and  has 
no  information  about  the  spectrum  of  eigenvalues.  Thus, 
with  no  experimental  data,  one  must  calculate  every¬ 
thing  a  priori . 

One  would  first  of  all  be  able  to  obtain  E  =  Bu , 
that  is  the  upper  bound  to  the  ground  state,  by  means 
of  the  familiar  variational  principle.  Now  for  the  first 
excited  state  one  would  use  Bu  =  E  as  E  in  Eqs.  (50) 
and  (48)  and  thus  be  able  to  obtain  a  B^,  which  represents 
an  upper  bound  to  the  first  excited  state.  To  determine 
the  upper  bound  to  the  second  state  one  would  now  use 
B^1  as  E^  and  proceed  as  before  to  obtain  B ^  .  In  order 
to  obtain  a  lower  bound  to  ,  one  must  first  of  all 
calculate  an  upper  bound  to  B^,  in  the  same  manner 

as  for  B^  and  B^.  Now  with  B^  as  E^  Eqs.  (50)  and  (48) 
may  be  used  to  obtain  a  lower  bound,  B^,  to  .  Care 
must  be  exercised  in  this  last  step  however,  for  B^  is 
of  necessity  above  ,  thus  two  true  eigenvalues  will 
be  found  between  B^  and  B^.  This,  however,  violates 
the  conditions  necessary  for  Kato's  theorem  to  be  true. 
Thus  one  would  have  to  use  a  value  slightly  less  than 
B^  in  order  to  avoid  this  difficulty.  The  inequality, 

A  <  (E  -  E^  ^) "  E) ,  should  aid  greatly  in  detecting 
problems  like  the  one  just  mentioned.  For  since  B^ 
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£ 

is  above  W^,  the  calculation  should  try  to  find  a  B^ 
which  would  give  an  E  above  .  But  in  this  case  when 
the  proper  quantities  are  substituted  into  the  inequality 
a  positive  quantity  on  the  left  side  and  a  negative 
quantity  on  the  right  side  would  result.  Hence  one  would 
know  that  something  was  wrong  and  thus  decrease  the 
value  of  B^.  This  is  a  very  simple  scheme  presented 
here;  more  elaborate  and  accurate  ones  certainly  could 
be  imagined.  Perhaps  in  the  future  it  would  be  worth 
while  to  investigate  in  detail  such  possible  alternatives, 
and  thus  be  freed  from  using  experimental  data  in  the 
calculations . 

I . 5  Considerations  of  Previous  Work 

There  have  been  few  applications  of  either  Weinstein’s 
or  Temple’s  formula  in  the  literature,  mostly  due  to 
the  difficulty  of  evaluating  the  integrals  in  A.  The 
few  applications  which  have  appeared  (7,  9,  10,  11), 
have  been  mainly  to  the  ground  states  of  He  and  H* . 

In  these  cases  the  $  which  is  used  to  calculate  the  lower 
bound  is  calculated  by  means  of  the  variational  principle, 
i.e.,  minimization  of  E;  in  Goodisman’s  work  on  H* ,  the 
minimization  of  a  is  used.  In  the  cases  where  E  is 
minimized  the  lower  bounds  are  poor,  but  according  to 
the  analysis  of  the  last  section  this  is  to  be  expected. 

For  if  $  is  obtained  with  the  criterion  of  the  minimum 
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E,  then  A  will  increase  and  since  A  or  its  square  root 
appear  in  the  numerator  of  both  formulae  this  will 
increase  the  quantity  which  is  added  to  E  in  both  cases. 
This  will  make  the  bound  worse.  In  the  case  of  A  minimiza¬ 
tion,  E  increases  and  thus  the  lower  bound  should  be 
better  than  in  the  case  of  E  minimization.  However 
neither  minimization  will  give  a  better  bound  than  the 
procedure  for  best  bounds  which  was  described  in  the 
last  section.  An  illustration  of  this  will  be  given 
in  Chapter  III  for  the  case  of  H*  . 

There  seems  to  be  some  confusion  in  the  literature 
on  bounds.  For  instance,  in  a  paper  by  Froman  and  Hall 
(7) ,  one  finds  the  statement  that  for  a  particular 
approximate  function  for  the  ground  state  of  He  the 
Weinstein  formula  gives  a  better  lower  bound  than  the 
Temple  formula.  This,  of  course,  is  incorrect;  the 
Weinstein  formula  can  never  give  a  better  bound  than 
the  Temple  equation,  as  was  shown  by  Kato  (3).  If 
one  examines  Table  I  in  their  paper  it  is  readily  apparent 
that  their  values  of  E  -  /a  and  E  +  /a  bracket  two 
eigenvalues.  This  contradicts  Weinstein's  initial 
assumption,  Eq.  (5),  and  hence  FrOman  and  Hall's  observa¬ 
tion  about  is  entirely  invalid. 

Goodisman  (12)  claims  to  have  derived  a  new  method 
for  obtaining  upper  bounds  to  excited  states.  Starting 
from  a  consideration  of  Stevenson's  formula  (13)  for 
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lower  bounds,  which  is  poorer  than  Temple's  equation, 
he  obtains  a  formula  for  upper  bounds,  which,  however, 
is  a  special  case  of  the  Temple-Kato  formula,  Eq.  (15). 

It  cannot  give  as  good  a  bound  as  Eq .  (15),  (3). 

In  a  recent  publication,  Goodisman  (5)  considers 
the  minimization  of  E  and  A  as  criteria  for  obtaining 
a  good  approximation  to  the  ground  state  eigenfunction 
of  the  system.  He  considered  the  variational 
function  of  Kolos  and  Wolniewicz  (14)  to  represent  the 
exact  function  closely  and  used  James  and  Coolidge  type 
functions  of  varying  complexity  as  trial  functions  . 

Then  after  obtaining  the  parameters  in  the  4^  by  either 
minimization  of  E  or  minimization  of  A  he  compared  the 
resulting  to  the  Kolos  and  Wolniewicz  function  by 
calculating  the  overlap  between  the  two  functions.  On 
the  basis  of  maximum  overlap  he  then  was  able  to  come 
to  conclusions  about  the  goodness  of  the  approximate 
functions.  His  results  showed  that  minimizing  A  did 
not  give  a  larger  overlap  with  the  "true"  function 
than  minimizing  E. 

Kolos  and  Roothaan  (15)  and  Kolos  and  Wolniewicz  (14) 
also  calculated  a  number  of  expectation  values  of  properties 
with  their  function  for  H^ .  Goodisman  chose  to  calculate 

<H>  ,  <3z2  -  r2>,  <l/r22>  an<^  <<S^12>  ^or  trial  functions, 

which  were  obtained  by  both  minimization  of  E  and  minimiza¬ 
tion  of  A  and  compare  them  to  the  "true"  values  obtained 
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by  Kolos  et.  al .  <H>  is  the  total  electronic  energy, 

<3z2  -  r2>  is  the  electronic  contribution  to  the 
quadrupole  moment,  <l/r^2>  Is  the  interelectronic  repul¬ 
sion  energy  and  <6r^2>  represents  the  size  of  the  set 
of  points  where  two  electrons  are  coincident.  The 
latter  term  arises  in  the  relativistic  correction  to 
the  energy.  and  <<$r^2>  be  particularly 

sensitive  to  local  fluctuations  of  the  wave  function; 
hence  one  would  suspect  that  they  might  be  better  calculated 
with  a  $  obtained  from  A  minimization.  However, the 
conclusion  from  the  data  presented  is  that  the  expectation 
values  calculated  using  the  functions  obtained  by  the 
minimization  of  E  are  closer  to  the  correct  values  than 
those  obtained  by  minimization  of  A.  The  exception 
is  <(5r-^2>  which  is  certainly  a  "very  short  range" 
property  and  highly  sensitive  to  local  energy  errors. 

It  is  a  bit  discouraging  however, that  <l/r^2>>  which 
should  be  sensitive  to  local  energy  errors,  is  better 
determined  by  the  E  than  by  the  A  minimized  $.  Goodisman's 
results  here  definitely  support  the  analysis  given  in 
section  1.3,  and  support  the  contention  that  a  minimization 
does  not  seem  to  be  a  good  method  to  determine  approximate 
wave  functions. 

A  general  method  for  obtaining  a  within  the 
orbital  approximation,  by  minimization  of  A,  has  been 
proposed  by  Fraga  and  Birss  (16) ;  since  the  method  is 
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basecl  on  the  minimization  of  A,  it  has  the  theoretical 
disadvantages  pointed  out  already.  Their  actual  procedure 
is  very  useful  however  and  much  of  their  work  can  be 
used  to  formulate  a  more  general  procedure  based  on  the 
arguments  presented  in  section  1.4.  It  will  be  shown 
in  the  next  chapter  that, using  their  basic  approach  of 
determining  pseudo-eigenvalue  equations  by  means  of 
coupling  operators  from  a  general  variational  principle  and 
combining  this  with  Eq .  (54),  one  is  able  to  develop 
equations  which  give  a  unified  approach  to  both  the  best 
bounds  and  the  best  wave  function.  These  equations, 
which  will  determine  the  best  bounds  or  the  best  wave 
function,  will  be  identical  in  form  and  will  differ 
only  by  the  constant  oo  which  will  appear  in  the  general 
coupling  operators . 

All  of  the  work  discussed  so  far  and  indeed  also 
the  method  which  we  will  develop  in  the  next  chapter 
are  based  on  some  sort  of  variational  procedure.  There 
are  also  methods  based  on  perturbation  theory  for  the 
determination  of  lower  bounds  but  these  methods  are  not 
capable  of  giving  upper  bounds;  thus  the  upper  bounds 
must  be  determined  by  another  method.  These  methods 
will  be  briefly  discussed  in  the  next  section. 

1.6  Other  Methods  of  Determining  Lower  Bounds 


The  first  method  of  determining  lower  bounds  which 
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will  be  considered  is  based  on  the  so-called  method  of 
intermediate  Hamiltonians.  The  approach  of  Bazley  (17) 
is  followed  here. 

The  method  consists  of  writing  the  total  Hamiltonian 
as  a  sum  of  two  parts 

H  =  H°  +  H\  (55) 

The  first  part  must  have  a  simple  enough  form  so  that 
the  Schrodinger  equation  involving  it  may  be  solved 
exactly.  The  second  term,  H',  must  be  a  positive  definite 
operator,  i.e.,  <  j  H 1  |  <f>>  >  0  must  hold  for  all  <f>  in 
the  domain  of  H’.  The  eigenvalue  problem  for  H°  is 
known  as  the  base  problem  and,  if  W?  are  the  eigenvalues 
of  H°,  they  will  give  rough  lower  bounds  to  the  eigen¬ 
values  of  H.  That  is,  one  may  write 

W°  £  i  =  1,  2,  3,  .  (56) 

this  is  a  consequence  of  the  operator  H’  being  positive 

1c 

definite.  Intermediate  Hamiltonians,  H  ,  are  constructed 
such  that 

H°  <  Hk  <  Hk+t  £  H  k,  t  =  1,  2,  3,  .  (57) 

The  operator  inequalities  in  Eq.  (57)  mean  that  for 
operators  A  <  B,  <<J>|A|<f>>  .<  <  4>  I B 1 >  for  arbitrary  <J>  in 
the  domains  of  A  and  B.  Bazley  (17)  and  Bazley  and  Fox 
(18)  have  given  some  ways  that  one  may  use  to  construct 
intermediate  Hamiltonians.  Then  as  a  result  of  an 
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ordering  theorem  they  obtain 


i.rO  k  -  Ch+  t  t.T 

W-  <  E  .  <  E  .  <W« 

1—1—1  ~  X 


k,t  =  1,  2,  3,  ...  (58) 


1c  1c  "t 

where  the  E^  and  E^  are  eigenvalues  to  the  intermediate 
Hamiltonians.  Hence,  as  one  constructs  intermediate 
Hamiltonians  which  approach  the  true  Hamiltonian,  one 
will  obtain  E^  which  serve  as  lower  bounds  and  approach 
the  true  .  Two  things  should  be  noted  however: 
firstly,  the  intermediate  Hamiltonians  must  be  so  con¬ 
structed  that  their  Schrodinger  equations  can  be  solved 
exactly;  and,  secondly,  if  one  wishes  to  obtain  a  lower 

bound  to  the  m^  state,  one  must  obtain  the  bounds  of 

t  h 

all  states  up  to  the  in  state.  The  latter  is  a  con¬ 
sequence  of  the  ordering  theorem. 

An  alternative  to  the  intermediate  Hamiltonian 
approach  has  been  suggested  by  Ldwdin  (19) .  His  approach 
involves  a  so-called  bracketting  theorem  by  means  of 
which  a  lower  bound  can  be  determined  to  an  eigenvalue 
if  one  has  an  upper  bound.  He  derives  the  equation 


=  <4>|H  +  H[P/(£  -  H)  ]  H  |  <J» 


(59) 


where  $  is  a  trial  function,  H  is  the  Hamiltonian  and 
the  quantity  P/(£  -  H)  is  the  reduced  resolvent  of  the 
problem.  If  one  knows  an  upper  bound  t  to  the  eigen¬ 
value  and  uses  equation  (59)  one  then  obtains  a  lower 
bound  t , .  Actually  Eq.  (59)  is  written  in  a  symbolic 
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notation;  to  write  it  out  is  rather  complicated.  H 
is  first  written  as  H°  +  V  and  considerations  of  infinite 
order  perturbation  theory  are  needed,  along  with  many 
operator  inequalities.  A  brief  outline  of  the  method 
is  presented  in  Appendix  III  in  conjunction  with  a 
calculation  presented  in  the  next  chapter.  Suffice  it 
to  say  here  that  the  method  avoids  the  problem  of  the 
ordering  theorem,  and  hence,  in  theory  at  least,  can 
be  used  to  determine  a  lower  bound  to  the  mL  state 
directly . 

It  should  be  noted  here  that,  unlike  the  method  of 
intermediate  Hamiltonians,  the  determination  of  a  lower 
bound  depends  upon  the  knowledge  of  an  upper  bound. 

The  better  the  upper  bound  is  known  the  better  the  lower 
bound  will  be.  This  has  been  shown  by  Reid  (20).  Choi 
and  Smith  (21)  have  recently  applied  both  the  method 
of  intermediate  Hamiltonians  and  the  Lowdin  procedure  to 
the  ground  state  of  He,  using  the  same  basis  sets,  and 
found  that  the  former  yields  a  better  lower  bound. 
Applications  of  Lowdin' s  method  to  excited  states  have 
not  as  yet  been  published. 

It  is  rather  difficult  to  see  how  either  of  these 
methods,  especially  Lowdin's  method,  can  be  applied 
simply  to  more  complex  systems.  Methods  which  depend 
upon  inner  projections  of  operators  onto  finite  subspaces 
seem  to  suffer  from  a  lack  of  knowledge  of  the  proper 

functions  to  span  the  subspaces  in  all  but  the  simplest 
examples . 
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I X .  Mathematical  Development 

_  4 

II . 1  Formulation  in  Terms  of  Orbitals 

In  this  chapter  a  general  scheme  based  on  Eq.  (1.54) 
will  be  developed  in  the  context  of  orbital  theory.  It 
is  assumed  here  that  the  trial  function  $  is  an  anti¬ 
symmetrized  product  of  one  electron  orbitals  or  a  linear 
combination  of  such  products,  such  that  4  possesses  a 
definite  orbital  and  spin  angular  momentum.  The  latter 
means  that  $  must  be  an  eigenfunction  of  the  L  ,  and  Sz 
operators.  One  term  in  4  may  thus  be  represented  as 

A[*1Cl)«(l)+1(2D8(2)4i2(3)a(3)*2(4DB(4)*3(5)a(5)  ...] 

(1) 

where  A.  represents  the  normalized  antisymmetrizer ,  the 
<j>^  are  the  space  functions  and  the  a(i)  and  B(i)  are 
the  spin  functions.  The  numbers  in  parentheses  represent 
the  coordinates  of  the  electrons  in  space  and  spin  space. 

If  K  is  the  non-relativistic,  spin  and  time 
independent  Hamiltonian,  then  the  expectation  value  of 
with  respect  to  $  is  written  as 

E  =  <K>  =  <$|k|$>  =  <<HI  Hu  +  ^  l  l/ra7r  |  4>>  (2) 

it  ofir 

where  it  is  assumed  that  the  function  $  is  normalized 
so  that  <$|$>  =  1.  The  summations  in  Eq.  (2)  are  over 
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all  electrons . 

The  expectation  value  of  ft  may  be  written  in  the 
general  case,  in  terms  of  orbitals  as, 


E  =  2  E  f*H>^  +  £  f*fyN 
£  x  11  -  - 


->->  1  m  limm 

lm 


(3) 


where,  as  in  the  work  of  Birss  and  Fraga  (16),  f^  is 

t  h. 

the  fractional  occupancy  of  the  l  shell  of  symmetry 

2 

designation  x.  For  example,  the  2p  configuration 

l 


wi 


ill  have  f £  =  fP  =  1/3.  The  £ 


l 


means 


,  that  is, 


i  Xa 


summation  over  all  orbitals  of  symmetry  species  x  and 
subspecies  a  and  summation  over  all  such  symmetry  species 
and  subspecies.  The  forms  of  and  along  with 

a  more  thorough  discussion  of  Eq.  (3),  are  given  in 
Appendix  I . 

It  will  be  recalled  that 


a  =  <;k>)2|$>  =  -  E2|$>  ;  (4) 


thus  one  must  also  obtain  an  expression  for  <$|lt2|$>  in 
terms  of  orbitals. 

The  operator ^f2  is  given  by 

•Ji2  =  I  HnH  +  2  I  l  ip2-  +  l  l  r-  C5) 

TT  O  TT  P  <  X  px  TT  <  O'  P<X  TT  (J  PX 

Now  evaluating  <$|ft2l$>,  where  Ti 2  is  given  by  Eq.  (5) 
and  subtracting  from  it  the  square  of  Eq.  (3)  one 


arrives  at 
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A  -  2  Y  £^Xr  t  A?+  +  T  B 

L  x  1  “  ■  X  m  im  “  .  x  m  n 


1 


-j~> 

xm 


imn 


xmn 


ii; 

lmnp 


f*fyfvf? 
1  m  n  p 


C*?— >■->•-> 
lmnp 


(6) 


where  the  A->+,  and  are  linear  combinations 

im’  xmn  lmnp 

of  integrals  which  are  given  in  Appendix  II  along  with 
a  derivation  of  Eq.  (6). 

It  will  be  recalled  that  the  $  which  satisfies 


6Z  =  6E  +  o)6A  =  0 


(7) 


will  determine  the  best  bounds  if  Z  =  B  and 
(E  -  Ey) 


03  = 


03  = 


l(E  -  E  r  -  A 

E  -  wkr 


2A 


or  will  be  the  best  if  Z  =  Q  and 


In  the  orbital  approximation  this  means 


'that  the  orbitals  which  determine  <±>  must  be  obtained 
and  that  they  must  be  such  that  the  $  satisfies  Eq.  (7). 

The  determination  of  such  orbitals  will  now  be  investigated 
The  variation  in  the  energy  in  terms  of  orbitals 
can  be  obtained  from  Eq.  (3)  and  is  given  by 


6E  =  2  £  {  <  6  <J>:f  |  f  ^  [  2H  +  £  f  yN->]  |  <p+>  +  c .  c .  } 
l  Ti'iL  L  m  mJ  1  1 


(8) 


-*■ 

x 


-+■ 

m 


The  c.c.  represents  the  complex  conjugate  of  the  first 
term.  Likewise  the  variation  in  A  which  can  be  obtained 
from  Eq.  (6)  is  given  by 
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-> 
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6A  -  l  {  <  6  $•>  |  f  ^  1 2X  +  £  +  l 


m 


->-> 

jnn 


fyfvB->+ 
m  n  mn 


fP  £V  f  ^ 

m  n  p  mnp 


]  |  <pf>  +  c  .  c .  } 


(9) 


Combining  Eqs .  (8)  and  (9)  according  to  Eq.  (7)  the 
following  is  obtained, 


6Z  =  l  {  < 5  <f)-^  |  f ^  [ 2  (H  - 

1 


A  A  A 


m 


+  c  .  c  .  } 


0 


(10) 


where  the  carets  over  various  symbols  denote  operators. 

In  obtaining  equations  (3)  and  (6)  it  has  been  assumed 
that  the  orbitals  are  orthonormal;  thus  one  must  constrain 
the  variation  of  the  orbitals  in  Eq .  (10)  so  as  to  maintain 
the  orthonormality  between  them.  This  is  accomplished 
by  the  technique  of  using  Lagrangian  multipliers. 

Originally  the  orbitals  are  such  that 


6  . 
1U 


(ID 


but  during  a  variation  they  must  satisfy 


P  ,  Xa i  .  Xa  '  .Xai^.Xa  n 
<6  d> .  <J>  >  +  <<t.  6  d>  >  =  0 

1 yu  1  Tu 


(12) 


in  order  to  remain  orthonormal.  Notice  that  only 
orbitals  of  the  same  symmetry  need  be  considered,  because 
two  orbitals  of  different  symmetries  are  automatically 
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orthogonal.  If  the  Lagrangian  multiplier  for  the 

constraint  involving  orbitals  i  and  u  is  designated 
X  a 

as  -  0^  and  this  is  multiplied  by  Eq.  [12)  and  one 
takes  all  such  products  and  sums  then,  one  obtains 


0 


(13) 


If  one  designates  the  operator  in  Eq .  (10)  as  Fj,  then 
that  equation  may  be  rewritten  as 


6 Z  —  £{<6c{>j|Fj|  4>j>  +  < S  |  F-£  |  <(>•£>  }  =  0 


(14) 


-> 

l 


Now  adding  Eq .  (13)  to  Eq .  (14),  one  obtains 


{  <6  d)^- 1  Ft)  <fr^ 
a1  Yi 


-  l 


\  n  \  n  *  .  +  * 

AaQA«>  +  <6,^  pT-j 

U  U1  1 


U 


l  **Aae4“>}  =  o 

L  TU  1U 

u  (15) 


Equation  (15)  will  hold  for  any  infinitesimal  variation 
* 

and,  since  the  variations  among  various 

orbitals  are  arbitrary  and  independent,  every  term 

in  curly  brackets  in  the  sum  of  Eq .  (15)  must  be  zero; 

* 

furthermore,  since  54^-  and  5  <}>-+  are  independent  of  one 
another  for  a  given  orbital,  each  of  the  two  integrals 
in  curly  brackets  must  independently  be  zero.  The  only 
way  the  integrals  can  be  zero  is  if  either  the  bra  or 
the  ket  is  zero;  the  bra  cannot  be  zero  since  arbitrary 
variations  have  been  assumed;  hence  the  ket  must  be  zero 
and  one  may  write 


■ 


5' 
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1  * 

u 

l  * 

u 

Ac(q  Xi£ 

u  ux 


* Aan Aa 
0-> 

U  1U 


(16) 


(17) 


Since  the  Lagrangian  multipliers  may  be  considered  as 
the  elements  of  a  Hermitian  matrix,  Eqs .  (16)  and  (17) 
are  equivalent.  If  the  system  under  consideration  is  a 
closed  shell  system  then  the  operator  designated  as  F-* 
is  identical  for  all  orbitals  and,  hence,  need  not  be 
subscripted.  Thus  one  may  write 


=  i  yv? . 

Yi  L  Yu  ui 
u 


(18) 


If  the  orbitals  are  subjected  to  a  unitary  transformation 
so  as  to  diagonalize  the  matrix  of  Lagrangian  multipliers 
the  following  equation  is  obtained,  since  the  operator 
F  will  remain  invariant  under  such  a  transformation 


Fd)*>  —  d)->0-}--^ 

1  i  11  ' 


(19) 


This  equation  is  well  known  (22)  for  the  case  that 


F  =  2H  +  y  N-> 

^  m 

m 


(20) 


Here,  however 

F  =  2  (H  -  toX)  +  y  (N->  +  wA:i) 

m 

+  l  u^mn  +  l  ^mnp)  C21) 

-y-y  -y 

mn  p 


. 

■ 
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for  the  case  of  closed  shells.  Thus  with  the  operator 
given  by  Eq.  C21) ,  one  finds  the  solutions  of  Eqs .  [19) , 
and  the  solutions  of  these  equations  will  be  such  that 
the  resulting  4  has  the  properties  previously  discussed. 
If,  however,  the  system  under  consideration  is  an  open 
shell  system,  then  one  must  follow  the  method  of 
Birss  and  Fraga  (16)  and  define  a  general  coupling 
operator  R 

This  operator  is  given  by 


where 


R 


Aa 


-  I 


Aa 


(22) 


<^>{><(*)5 1 F£)  "  ^  Uj><<f>glF 

V 


Aa i . Aa  ,  A  a 

(D  >  <  m 

v  1  Y  V  ' Y  V 


<(f>£  I  Fg  I  <(f>£ 


- i 

v 


,  Aa  .  Aa 
<b  >  <  d) 

TV  yV 


,A  a 


Fv  I 


+ 


<(P 


(23) 


In  a  general  open  shell  case  the  operator  F£  which  appears 
in  Eq.  (23)  is  given  by 


Fy  =  f,  {2(H  -  wX)  +  T  fy(N+  +  U)A+)  +  y  fyfVco 
k  k  ^  J  imvm  nr  A  m  n 

m  mn 


x  (B++  + 
v  mn 


7  f^O++)}. 

*  p  mnp' 

P 


(24) 


By  use  of  Eq .  (23),  Eq .  (16)  may  be  written  in  the  form 


nAa  ,  Aa  ,  Aan Aa 

R  0^  • 


(25) 


.-*V  - 


. 
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Hence  the  Eqs .  (16)  have  been  transformed  into  pseudo¬ 
eigenvalue  equations  which  may  be  solved  by  standard 
iterative  techniques.  There  will  be  a  set  of  equations 
like  Eqs.  (25)  for  every  group  of  orbitals  of  symmetry 
designation  Xa,  y 3,  ...  .  After  having  determined 

all  these  groups  of  orbitals  from  the  solutions  of  Eqs. 
(25),  one  may  then  form  $  from  these  orbitals.  This 
$  will  be  the  desired  solution  of  Eq .  (7)  and  will  hence 
serve  to  give  the  best  bound  or  best  4^  depending  upon 
the  quantities  Z  and  to . 


1 1. 2  Expansion  of  the  <f>^  in  Terms  of  a  Basis  Set 


X  a 

The  individual  orbitals,  <|k  ,  may  be  expanded  in 

X  c t 

terms  of  a  given  basis  set  of  functions,  >  as 


,  X a  r  Xa^X 
<p  •  =  >  x  c  - 

l  L  q  qi 


(26) 


where  the  are  the  expansion  coefficients.  Eq .  (26) 

may  also  be  written  more  conveniently  in  matrix  notation 


as 


,Xa  -*Xa-±X 
♦i  =  x  ci 


(27) 


where  x"  is  a  row  vector  whose  elements  are  the  basis 

functions  of  symmetry  designation  Xa  and  is  a  column 

vector  whose  elements  are  the  coefficients  . .  The 

qi 

latter  are  independent  of  subspecies  designation,  as  they 
must  be  in  order  to  maintain  the  correct  symmetry 


* 


. 


'  - 


■ 


. 


■ 
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transformation  properties  among  the  orbitals.  Since  the 
orbitals  have  been  assumed  orthonormal  this  leads  to 


Cpnj0(p)dvp 


X  e  X  X  _  . 

•  w  Ij  •  0  «  • 

j  ij- 


(28) 


X  c t 

where  S  is  the  overlap  matrix.  The  elements  of  the 
overlap  matrix  are  given  by 


The 

the 


c  X  a  Xa  |  X  a  r  X  o t  ^  ^  X  a  ■v  i.t 

Sqr  <xq  I  *r  *  "  >  xq  (pUr  (p)dVp. 


eigenvalue  equations  (25)  may  now  be  expressed 
form 

0 


RXa2x  =  ex“sXa2^ 


ii. 


(29) 


in 


(30) 


where  RXa  is  the  matrix  representation  of  the  operator 
RXa,  i.e.,  RXa  =  <xXa|RXalxXa>>  in  terms  of  the  basis 
set  of  functions.  There  will  be  such  a  set  of  equations, 
Eqs .  (30),  for  every  symmetry  designation  Xa. 

In  the  case  of  a  $  which  has  M  occupied  orbitals 
of  type  Xa  and  N  occupied  orbitals  of  type  yg,  etc.  [for 
the  upper  bound  and  best  calculations]  the  first  M 
orbitals  of  type  Xa  and  the  first  N  orbitals  of  type  yg, 
etc.  will  be  taken  from  the  solutions  of  Eqs.  (30)  in 
order  to  obtain  the  proper  $.  This  is  because,  in  these 
two  cases,  one  is  minimizing  Z,  i.e.,  minimizing  Q2  in 
the  case  of  finding  the  best  and  minimizing  the  upper 
bound,  Bu,  in  the  other  case.  However,  in  finding  the 
lower  bound,  B  ,  one  is  maximizing  B  and  hence,  for  a 


0 


- 


■ 
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$  as  described  aboye,  one  would  take  the  last  M  orbitals 
of  type  Aa  and  the  last  N  orbitals  of  type  yB,  etc.  in 
order  to  obtain  the  appropriate  $.  This  problem  of 
picking  the  right  set  of  orbitals  is  of  course  a  consequence 
of  Eq.  (26)  ,  since  if  there  are  M  orbitals  of  type  Aa  the 
sum  over  q  may  run  up  to  q  =  U  where  U  >  M. 

In  terms  of  the  basis  functions,  the  integrals  of 
Appendix  II  have  the  following  forms 


1  =  /  X*(p)Hpxs(p)<lVp 
II  -  /  [Hp+x*(p)][HpXs(p)]dVp 


(31) 


(32) 


III  =  /  /  X*(p)xs(p)  d/rpt)x‘(T)xw(T)dVpdVT  (33) 

IV  =  /  /  X*(p)xr(p)(l/rpT)x*(T)[HTxt(T)]dVpdVT  (34) 

V  =  /  /  Xq(p)xr(p)(l/rpT)2x*(T)xt(OdVpdVT  (35) 

V!  =  /  /  /  xq(p)xr(p)(1/rpT)xs(T)xt(x)(l/rpn) 


x  x*(n)xvCn)dVpdVTdVri. 


(36) 


Here  the  subscripts  on  the  basis  functions  denote 


different  basis  functions  of  the  same  symmetry  o_r  of 
different  symmetries,  i.e.,  q  may  be  of  type  Aa  and  r 
of  type  Aa  or  yB,  etc. 


- 


- 
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The  correspondence  between  the  integrals  above  and 
the  ones  in  Appendix  II  is 

I  ->  H 
II  ->  X 
III  +  N 
IV  ->  P 
V  ->  Q 

VI  +  T  and  Y 

where  it  is  to  be  understood  that  linear  combinations 
of  the  integrals  on  the  left  give  the  type  of  integral 
denoted  on  the  right.  The  most  common  type  of  basis 
function  used  in  the  expansion,  Eq .  (26),  is  the  Slater- 
type  -orbitals  (STO’s).  A  discussion  of  the  various 
integrals  which  arise  from  the  forms  I  -  VI  for  one  and 
two  center  problems  using  STO’s  will  now  be  given. 

I I . 3  Integrals  over  Slater  Basis  Functions 

One  Center  Case 

For  the  one  center  case  the  one  electron  operator 

is  given  by  V2  z 

H  =  -  -p- -  (37) 

p  2  rp 

where  Z  is  the  nuclear  charge  and  r  is  the  distance  of 
til 

the  p  electron  from  the  nucleus.  The  two  electron 
operator  is  1/r  ,  where  r  is  the  distance  between 

r  pi  p  T 

electrons  p  and  x  . 


i 

. 


. 
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Evaluation  of  the  integrals  I  and  III  for  this  case 
is  quite  familiar  and  the  techniques  are  described  in 
standard  text  books  C23) .  The  evaluation  of  integrals 
II  and  IV  for  the  one  center  case  can  be  accomplished 
by  the  same  methods  used  to  evaluate  I  and  III  with 
certain  obvious  modifications;  the  result  of  IV  is  a 
linear  combination  of  integrals  of  type  III.  However, 
the  values  of  the  integrals  V  and  VI  for  the  one  center 
case  are  not  so  easily  obtained  and  the  method  of 
evaluation  not  so  obvious.  An  analytic  evaluation  of 
V  was  accomplished  by  Solony,  Lin  and  Birss  (24)  ,  and 
the  analytic  evaluation  of  VI  will  be  given  here. 

As  in  the  work  of  Solony,  Lin  and  Birss,  unnormalized 
Slater  functions  are  used  here: 


n-1  -5rn  m 
e  P' 


(cos  0)  f  (  |  m  |  <j>)  (38) 


(n,  i,  m)  =  r 


Si 


f  =  cos,  for  m  >_  0 ;  f  =  sin,  for  m  <  0.  The  integral 
in  question  may  be  rewritten  as 


<(n1,  <4  >  mP  1 1/rlc Knc >  lc’  mP 


-III  X1COxiCl)Cl/rlc)xcCc)x'CcKl/r2c) 


.  / 


X2C2)x|  C2)dV1dVcdV2 


C39) 
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where  the  subscript  and  electron  label  c  emphasize  the 

role  of  these  functions,  i.e.,  electrons  1  and  2 

interact  with  electron  c  but  not  with  each  other. 

The  products  Xl  Cl)  Cl)  ,  x2C2)x2C2)>  xc  (c)  X(L  (c)  , 

can  be  separately  expanded  as  linear  combinations  of 

unnormalized  charge  distributions  [N-^ ,  ,M^  ]  ,  [N^  ,  ,1^  ]  > 

[N  ,L  ,M  ]  with  the  result  that  the  integral  VI  becomes 
c  c  c 

a  linear  combination  of  integrals  over  basic  charge 
distributions,  of  the  form 


(40) 


with  l±  =  (e ±  +  x,  j)/2  . 

Two  reasonably  distinct  methods  of  evaluating  these 
integrals  have  been  investigated.  In  one  of  these,  the 
integrals  over  electrons  1  and  2  in  Equation  (40)  are 
evaluated  by  the  scheme  giyen  by  Roothaan  (25),  yielding 
functions  of  electron  c.  Combination  of  these  functions 
with  the  charge  distribution  for  electron  c,  and  evaluating 
the  remainder  of  Equation  (40)  leads  to,  among  others, 


■ 
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the  integral 


/ 


Nc-Li-L2-1  _-2C€c  +  4  +  C2]rc 


0  rc 


dr 


Unfortunately  there  are  cases  of  interest  for  which 
N  -Ln-L?-1  is  negative.  In  the  second  method,  1/r.  is 
expanded  in  spherical  harmonics  and  the  radial  variables 
r  and  r  ,  the  lesser  and  greater  of  r-  and  r  .  This 
requires  that  the  integration  limits  for  r^  and 
depend  upon  r  .  Comparison  of  the  two  methods  led  to 
the  conclusion  that  the  difficulty  essentially  arises 
because  one  introduces  an  infinity  when  evaluating, 
e.g.  , 

/0C  f(r1)dr1  =  /q  f(r1)dr1  -  /”  fCr1)dr1  (41) 

c 

a  technique  which  leads  to  forms  identical  to  those  of 
Roothaan  at  the  stage  just  before  integration  over 
electron  c.  Therefore  a  modification  of  the  second 
technique  was  adopted;  it  avoids  the  introduction  of  the 
infinity  by  a  reversal  of  the  order  of  integration. 


Integration 


Th 


e  operator  l/r^c  is  expanded  as 


1/r. 


1C 


Jo  Jo  Mf 


r< 


r> 


k+1 


P^Ccos  ei) 


i,c 


P^Ccos  ec)  cos  m (4^  -  4c) 


(42) 


X 


- 
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where  r  and  r>  are  the  lesser  and  greater,  respectively, 
of  r-  and  r  .  Substitution  of  this  form  for  l/r1r 
into  Equation  (40)  introduces  the  factor 


|M.  I 

/ 0  PL1  (cos  V 


P^Ccos  e1)sin  e1  d0x  /^f  (  | M1 1  cf^) 


x  cos  -  <f>c)d(})1 

2,(!  +  6  )  (Ll  +  IMJ)! 

= _ _ _ g 

2L1  +1  (Lx  -  |M1 | ) l  kLl 

Only  the  term  with  k  =  ,  m  =  | 

and  that  with  k  =  L^,  M  =  |M2|  in 
(40)  becomes 


6m|M1 1  f ^ I Mi I ^  c  ^ * 

|  subsists  in  1/r^  , 
l/r2c,so  that  Equation 


U  -  16^2[9]U]  r  r  1 

U  (2L1+1) (2L2+1)  J  rl 


N.,+1  -2^2T1 

e 


r< 


Ll  +  ! 


c  r> 


dr- 


N  +1  -2^  r 

r  c  c  c 

X  /  rc  e 


,  r  N2  +  1  ~ 2  ^  2r2 

drc  /  r2  e 


"l,c 
r  L„ 


r  < 


V1 

r> 


dr. 


'2  c 

^  i 


(43) 


with 


|M, 


I M- 


!  M. 


[e]  '  /o  PL  1  (-C0S  6p  PL  2  (-COS  e2-)  PL  ^  (cos  ec) 
1  2  “  c 


x  sin  0  d0_ 

c  c 


(44) 


[4]  *  !2*  f(|M1lfc)f(|M2Uc)f(|McUc)  d*c. 


(45) 
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Four  distinct  cases  can  be  distinguished,  arising  from 
various  assignments  of  r  ,  and  r  : 


A: 

H 

A 

t-1 

V# 

O 

II 

H 

n 

v# 

r>l,c  = 

rn 

r<  2 , c  V 

r>  2 

B: 

r<l,c  =  rc> 

r>l,c  = 

ri; 

r<2 , c  =  r2 ’ 

r>  2 

C: 

r<l,c  =  rl’ 

r>l,c  = 

v 

r<2 , c  =  rc’ 

r>  2 

D: 

r<l,c  =  rl’ 

r>l,c  = 

v 

r<2 , c  =  r2 ’ 

r>  2 

c  "  r2 


=  r 
c  c 


=  r. 


=  r 

c  c 


For  case  A  the  radial  integrals  of  Equation  (43)  give 


N  +Ln+L0  +  1  -2 e  r  N-.-L-, 

I,  =  /“  r  c  1  2  e  c  c  dr  /”  i,1  1  e  1  1  dr. 

a  J  0  c  cJ  r  1 

c 


N0-L0  -2<;0r0 

r°°  2  2  2  2  , 

x  I  r0  e  dr, 

J  r  2 
c 


0). 


N..-L-,  ,  N0-L0 

'  1  1  k  2  2 

N,  - L,  + 1  N,-L,+l  N  +  N.,+N,  +  4  A-  FT  Jn 
11  2  2  c  1  2  Jc=0  £  =  0 


(N1-L1) ! (N2-L2) 


GO  . 


(Nc+Li+L2+k+£+l) !  £ 

-  0) 

A!  Z 


with 


For  case  B, 


=  2^/a,  a  =  2(cc  +  ^  +  ?2) 


(46) 


T  ..  Nc+LrLz  '2?crc  ,  r 

XB  =  / 0  rc  e  drc  It  rl 

-  c 


NrLi  -2?iri 
e 


dri/i 


N2+L2+1  '2c2r2 
x  r2  e  dr 2 


. 
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NrLi 


=  (N1-L1) ! (N2  +  L2  +  l) !  / 


l 


N.  -L  +1  N  ,  +  L ,, + 2  N" +'n“-'n . t'4 

11  2  2  cl2 

M11  11  21  P1  k=0 


CNc+LrL2+k) !  k 

-  y 


k! 


11 


N2‘L2  W1 


N,  -L-.  +  1  N0  +  L0  +  2  N  +N-.  +N0  +  4  ,  Ln  k! 

11  2  2  c  1  2  k=0 

a)  to  a 


W1 

I 

£=0 


(Nc+L2-L2+k+i)! 
£  ! 


0}. 


£  V 
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with  y .  .  =  2  £./p.,  p.  =  2  (e  +  ?.)• 
1J  1  1  ’  1  ^  c  1 


For  case  C  the  form  is  that  of  case  B  with  labels  1  and 
2  interchanged  throughout.  In  both  of  these  cases  use 
of  Equation  (41)  has  been  made  since  the  critical 
exponents  of  r  ,  N  h-L^-I^  or  N^^-L^  +  L^,  cannot  become 
negative  (this  is  a  result  of  the  inequality  given  in 
Equation  (5  2))  . 

In  case  D, 


N  -L-,  -L0-l  -2 1  r 

T  r00  cl  2  ^  c  c  ,  .  r 

XD  =  '0  rc  e  drJ 


r  NpL  +1  -Zcr 
q  r1  e  dr- 


rc  N2tL2+l  -2c2r2 

*  f 0  r2  6  dr2> 


this  is  the  only  term  giving  rise  to  the  difficulty 
discussed  in  the  introduction.  To  obviate  the  trouble 


. 

-Sl¬ 


one  can  change  the  integration  variables  and  then  the 
order  of  integration  to  obtain 


N  +N,  +N0  +  3  -2 1  r  N-.+L..+1  -2c, r  x 

T  r°°  C  1  2  C  C  ,  fl  Y  1  1  1  C  J  r  1 

ID  =  JO  rc  e  drc  Jo  x  e  dx  fo 

N2+L2+1  -2«2rc7  , 

x  y  e  dy 


,  N, +L, + 1  -|  N.  +  L.  +  l  N  +N,  +N0  +  3 

rl  x  1  1  j  r  1  22  ,  r00  c  1  2 

=  /  0  X  dx  Jo  y  dy  /  o  rc 


x  e 


-2(?c+?lX+'2y)r, 


dr 


(Nc+N1+N2+3) 


Nl+Ll+1 


N1  +  L1  +  l\ 


N  +N-.+N.  +  4  N  +Ln  +  2  N.  +  L.  +  2  ,  Ln 

?cl2  -11-22  k=0 

Z  ^1  ?2 


(-D 


N  -L1+N0+2+k 
c  1  2 


N?+L?+l/  v 

x  Y  (N2+L2+1 

(-if- 

-  J l 

II 

o 

*  L 

l/l 


N  -L, -L0+A 
c  1  2 


1/ (?c+^2^ 


N  -L-.-L-  +  JI 
c  1  2 


/(Nc"LrL2+£) 


m=0 


m 


N  -L, -L~+£+m 
c  1  2 


-i/Ccc+51+c2D 


N  -L, -L0+£+m 
c  1  2 


l 


/cVVl2+£+iu)**j  C48) 


The  asterisk  indicates  that  when  Nc-L^-L2+£  =  0,  it  and 
the  preceding  term  in  square  brackets  are  replaced  by 
In  [  Cc c+c /c  ] ;  similarly  for  the  double  asterisk,  if 
Nc-L-^-L2  +  ii+m  =  0,  it  and  the  preceding  term  in  square 
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brackets  are  replaced  by  In  [  (A  +£  S  2^  /  (4C  +  ^)]  • 

If  N  ,  and  are  such  that  N^^-L^-L^-l  is  positive, 
the  method  of  Equation  (41)  can  be  used  to  evaluate  1^, 


D 


CN1  +  L1  +  1)!CN2  +  L2  +  1.).!  1 

[Cfr-q-q-D! 

N  +Nn+fOT  N7+L7T2  N,+L,+2l 

2  c  1  2  q  1  1  q  2  2  1 

1.  VLrT2 

VV1  (Nc-Ll-L2+k-l)! 

x  1 
k=0 


N  -L-,  -L, 


11 


k!  ^c+Cl} 


N2+L2+1  (N  -L1 -L0+k-l) !  , 

Y  c  1  2  k 

1  n  N  -L-.  -L„y2 

k=0  k l(J  +52)  c  1  2 


N  +  L  + 1  k  Nn  +  L~  +  1  r .  T  T  T  , 

11  w.  2  2  (N  -L-. -L0+k+£-lJ  ! 

V  1  v  v  c  1  2  J  l 

Jo  - n - “2 


JL 
k! 


I 

Z  =  0 


(49) 


The  formulae  given  above  were  checked  in  a  number 
of  ways.  Using  the  alternate  form  for  1^,  Equation 
(49) ,  the  analytic  form  of  U  derived  here  was  checked 
against  a  number  of  specific  formulae  obtained  by  the 
use  of  the  Roothaan  scheme.  Two  versions  of  a  computer 
program  were  made,  using  Equations  (48)  and  (49)  for  1^. 
Where  (49)  was  applicable,  the  two  methods  agreed. 
Numerical  differentiation  was  used  in  conjunction  with 
the  program  using  Equation  (48)  to  raise  the  quantum 
numbers  of  the  various  charge  distributions.  Agreement 
between  values  obtained  in  this  way  and  by  calculation 
using  the  raised  quantum  numbers  in  Equations  (46)  to 
(48)  was  found.  In  particular,  in  those  cases  where  Nc 
was  raised  from  a  condition  where  the  logarithmic  terms 


. 
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entered  to  a  high  enough  yalue  that  they  no  longer  occurred, 
the  agreement  can  be  taken  as  guaranteeing  the  validity 
of  Equation  (.48)  . 

The  Angular  Integrations 


The  integrals  [e]  and  [4]  of  Equations  (44)  and 
(45)  are  well  known  in  the  theory  of  vector  coupling. 
[4]  is  zero  if  one  of  M-^ ,  M2,  Mc  is  negative,  or  if 
all  three  are  negative  or  if  the  equality 


M.  +  M.  =  M, 
1  1  1  3  1  1  k 


is  not  satisfied  by  some  assignment  of  1,  2,  and  c  to 
i ,  j  ,  and  k .  Then 


[  4>  ]  =  2 7T ,  M.  =  M.  =  Mk  =  0 

3 


=  TT,  M.  =  M.  +  0,  M,  =  0 

*  1  j  ’  k 


=  -tt/2,  M.  f  M.  ,  M-  <  0,  M.  <  0,  Mt  >  0, 
1  J  1  9  3  k  ’ 

M.  +  M.  =  M, 

1  j  k 

=  tt/2,  for  all  ?  0 ,  the  not 
satisfying  the  above  conditions. 


One  possible  form  of  [4]  is  C26) 


/  P*  (cos  4)  (cos  4)  P^  (cos  4)  sin  4d4 


< 
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=  (;  _  i)  C ^ -m+n)  / 2  - k 

x  _ (m+j)  !  (n+lc)  !(£+m-n)  !  C(£+m+n)/2)  I _ 

(m- j  )  !  (£+m+n+l)  !  ((£+m-n)/2)  !  (  (£-m+n)  /  2)  !  (  (-£+m+n)  /  2)  ! 

x  £  _ C^  +  i+p)  I  (m+n-i-p)  ! 

p  p! (fc-i-p) ! (n-k-p) ! (m-n+i+p) ! 

with  the  restrictions 

i+j=k,  m+n>£>m-n,  m>n.  (52) 

This  completes  the  discussion  of  the  evaluation  of 
integral  VI  and  with  its  evaluation  all  the  integrals 
necessary  to  carry  out  a  calculation  involving  A  on  any 
atom  are  now  available.  Hence  the  methods  described 
in  Chapter  I  and  developed  in  this  chapter  in  terms  of 
orbital  theory  can  be  applied  to  any  atom,  so  long  as 
the  orbitals  are  given  by  linear  combinations  of  STO's. 
However,  this  situation  does  not  exist  for  diatomic 
molecules,  as  will  be  shown  below. 


Two  Center  Case 


For  the  two  center  problem,  the  one  electron  operator 
is  given  by 

V-  ZA  ZR  *1  ZA  ZR 

E  B  A  B  (53) 


h  =-■-£■-  -A 


'A 


'B 


2  r  a  r  n 
p  A  pB 


2  rA  rB 


where  and  Zg  are  the  nuclear  charges  of  nuclei  A  and 
B  respectively,  and  where  r  ^  = 


rA  and  rPB 


rB  are 


- 
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the  distances  between  electron  P  and  nuclei  A  and  B 
respectively . 

Evaluation  of  integrals  I  and  II  for  the  two  center 
case  has  been  accomplished  [25,  27).  In  the  integral 
forms  I,  II  and  IV,  H  appears  where  x  is  a  STO  with 
orbital  exponent  a  and  H  is  as  given  in  Eq .  (53). 


Hx  = 


-n(n-l)  -  £(£+1)  + 

2r2 


a  - 


Z 


r 


(54) 


which  shows  that  the  one  electron  operator  acting  upon 
an  STO  produces  a  linear  combination  of  STO’s.  Thus  the 
integrals  having  the  form  I  can  be  expressed  in  terms 
of  the  familiar  overlap  type.  The  integrals  having  the 
form  III  are  of  three  types: 


(1) 

<XAXA 1 1/r12 

lxBxBf>  - 

coulomb  type 

(2) 

<XAXB  1  1//r12 

1 XBXA ' >  ’ 

exchange  type,  and 

(3) 

<xAxaI 1/r12 

1 XAXB ' >  ’ 

hybrid  type. 

Where  the  A  and  B  denote  the  nuclei  on  which  the  STO 
is  centered. 

In  considering  integrals  of  the  form  II,  the 
following  types  occur,  as  can  be  seen  by  using  Eqs . 
(53)  and  (54) , 


1 
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CD 


(2) 


(56) 


where  m,  n  =  0,  1,  2  and  m  +  n  <  2.  These  integrals 
have  been  evaluated  analytically  (24,  28). 

For  the  integrals  of  the  form  IV,  three  types  exist 

(1)  <xaxAU/  r12 I xBH^B ' >  *  a  coulomb  tyPe 

(2)  <XaXb I !/r12 I XgHxX' >>  an  exchange  type,  and  (57) 

(3)  <xAxkW*12h'PxB'>’  a  hybrid  bype- 

Each  of  these  integrals,  as  a  result  of  Eq .  (54), 
will  be  a  linear  combination  of  integrals  of  the 
corresponding  type  given  in  (55)  .  Hence  the  evaluation 
of  the  basic  integrals  involved  in  forms  I-IV  has  been 
accomplished.  However  the  forms  V  and  VI  are  considerably 
more  difficult  and  no  analytic  expressions  have  been 
derived  to  evaluate  these  integrals  as  yet.  The  various 
types  of  integrals  which  must  be  evaluated  as  a  result 
of  the  forms  of  V  and  VI  will  be  briefly  considered. 

For  the  integrals  of  form  V,  again  three  types  of 
integrals  will  arise: 


(1) 


* 


. 

■ 
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(2) 

<XAXB>  Cl/rlzJ2  |-xgxS'> 

(58) 

(3) 

<XAXaI  *-1/r12-*  2  1  xAXb’>  ■ 

The  difficulty  in  evaluating  these  integrals  arises  from 
the  (l/r-^)2  expansion  onto  two  centers. 

For  the  integrals  of  the  form  VI,  the  following 
types  arise, 


f'i'i  Iii/  |  II  III  i -i  /  |  IV  V 

C1)  <<^AXA  1 1// r12  I  XA  XA  l1/r13  I  XB  > 

C2)  <XAXA  l1/r12  I  XA  XB  1 1  r  1 3  I  XA  XA  > 

(3)  <XAXA  l1/r12  I  XA  XA  1 1  ^ r  1 3  I  XB  XB  > 

(4)  <XAXB  1 1^r12  I  XA  XA  1 1  ^ 17 1 3  ^  XB  Xa> 

rc>  I  M  ,  ,  II  III  ,  ,  IV  V 

C5)  <XAXA I 1^r12 I XB  XB  I 1^r13 I XA  XB > 

(6)  <XAXB 1 1 ^ r 1 2 1 XB  XB  1 1 ^ r 1 3  I x A  Xa> 

C7)  <XAXA 1 1 ^ r 1 2  I XA  XB  l1^r13lxB  XA> 

(8)  <xaXb  1 1//r12  I  Xa^  B  l1(/r13lxA  XB> -  . 


' 
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Thus  in  order  to  apply  the  methods  discussed  to  a 
diatomic  species  with  three  or  more  electrons  the  eleyen 
integrals  of  (58)  and  (59)  must  still  be  evaluated. 

This  is  obviously  a  place  where  further  work  is  desirable. 

Considering  the  difficulty  involved  in  evaluating 
three  and  four  center  integrals  of  H  over  STO's  the 
number  and  difficulty  of  integrals  involving  Ti 2  would 
seem  to  preclude  any  thought  of  extending  methods 
involving  H 2  to  polyatomic  systems,  until  a  thorough 
investigation  of  these  methods  on  one  and  two  center 
problems  has  been  carried  out. 

Thus  the  methods  described  in  the  preceeding 
sections  are  at  present  limited  to  atoms  and  one  electron 
diatomic  species . 
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III.  Calculations  on  the  H*  and  He  Systems 

1 1 1 . 1  Introduction 

In  section  2  the  ground  state  of  the  H*  system  is 
considered  and  results  are  obtained  and  compared  for 
the  lower  bounds  obtained  by  three  variational  schemes: 

E  minimization,  the  A  minimization  scheme  of  Fraga 
and  Birss  (16)  and  the  lower  bound  maximization  scheme 
described  in  the  last  chapter. 

In  section  3  upper  and  lower  bounds  are  obtained 
3  1 

for  some  ’  P(ls,np)  excited  states  of  He.  The  upper 
bound  results  are  compared  with  upper  bounds  obtained 
by  Davidson  (30)  for  these  excited  states. 

The  method  proposed  for  obtaining  the  best  <5^  for 

3  1 

an  excited  state  is  applied  to  the  *  P(ls,np)  excited 
states  of  He  in  section  4.  These  wave  functions  are 
tabulated  there. 

The  wave  functions  obtained  in  sections  3  and  4  are 
used  in  section  5  to  calculate  some  physical  properties 
of  the  He  excited  states.  These  properties  are  tabulated 
and  a  discussion  of  the  results  is  included. 

In  section  6  a  Hylleraas  type  calculation  is  carried 
out  for  the  ground  state  of  He  in  order  to  provide  a 
direct  comparison  between  the  method  proposed  by  Lowdin 
for  lower  bounds  and  the  method  proposed  in  this  work. 
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II I. 2  Calculations  on  the  Ground  State  . 


In  this  section  the  results  of  the  application  of 

>  *4* 

the  method  developed  in  Chapter  I  to  the  ground 
state  will  be  presented.  In  applying  the  formalism  for 
lower  bounds  the  most  practical  method  is  the  expansion 
of  the  orbitals  in  some  suitable  basis  set,  as 

mentioned  in  section  II.  2.  One  may  write 


Equation  (2.25)  c 

coefficients,  c. 

’  IP 

RXa 

operator  unde 
The  procedure  is 
self-consistency 
pseudo -eigenvalue 
the  operator  (or 
self-consistency 
are  then  used  to 
by  Eq.  (1.48)  . 


,  X a  =  r  X  Xa 
4>i  L  cii 


-P  XP 


CD 


an  be  expressed  in  terms  of  the  linear 
,  and  the  matrix  representations  of  the 
r  the  basis  set,  as  given  in  Eq .  (2.30) 
then  the  familiar  one  of  any  iterative 
approach,  with  the  solutions  of  the 
equations  being  used  to  construct 
its  matrix  representation) ,  until 
is  attained.  The  resultant  orbitals 
calculate  B,  the  lower  bound,  as  given 


It  should  be  noted  that  the  upper  bound  for  the 
ground  state  as  given  by  Eq.  Cl. 48)  reduces  to  the 
ordinary  variational  principle  when  E^  -  -°°,  which  is 
the  proper  choice  of  E^  for  calculations  that  are  to 
determine  the  upper  bound  for  the  ground  state.  Thus, 
since  the  upper  bound  in  this  case  is  determined  by 


I 

i\'  [f  ft  .]  3J  if  3  ; 


■ 


' 


. 
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the  variational  method  for  lowest  states,  only  the  lower 
bound  need  be  considered  here. 

The  trial  function  used  is  a  generalization  of  the 
James  function  (31) , 

*  =  e'a5  l  c.n2i  (2) 

•  _L 

l 

where  £  and  r\  are  the  variables  in  the  confocal  elliptic 

coordinate  system.  The  integrals  were  evaluated  using, 

in  part,  the  analytic  forms  of  Goodisman  and  Secrest  (9). 

Besides  the  method  developed  in  this  work,  functions 

obtained  by  the  scheme  of  Fraga  and  Birss  and  by  the 

ordinary  variational  principle  were  also  used  to  obtain 

a  lower  bound  to  the  ground  state.  All  three  methods 

are  used  to  obtain  the  molecular  orbital  <J>  for  the 

ground  state  and  these  (p's  are  then  used  in  the  Weinstein 

formula,  Eq.  (1.5)  and  the  Temple  formula,  Eq .  (1.48) 

in  order  to  obtain  the  lower  bound.  The  three  methods 

each  yield  a  different  <f>  and  therefore  give  different 

lower  bounds.  The  results  are  tabulated  in  Table  1. 

The  variational  schemes  designated  as  E,  A  and  B  are  the 

ordinary  variational  principle,  the  A  minimization  scheme 

of  Fraga  and  Birss  and  the  lower  bound  maximization  procedure 

proposed  in  this  thesis,  respectively.  It  is  clear 

from  the  results  that  the  scheme  proposed  here,  used 

in  conjunction  with  the  Temple  formula  gives  the  best 

lower  bound.  BT  ,  and  B1>r  .  .  .  give  the  bound  as 

Temple  Weinstein  a 


Table  1 


Comparison  of  Lower  Bounds 


No.  Basis 
Functions 

A 

g 

Weinstein 

^Temple 

a 

Variational 

Scheme 

1 

0.091676 

-1.38116 

-1 . 20615 

1.400 

A 

2 

0.001422 

-1.14009 

-1.10430 

1.355 

E 

2 

0.001272 

-1.13797 

(-1.14724)* 

-1.10403 

(-1.10706)* 

1.340 

A 

2 

0.001444 

-1.13994 

-1 . 10388 

1.320 

B 

3 

0.000903 

-1.13251 

-1.10367 

1.355 

E 

3 

0.000519 

-1.12511 

-1.10303 

1.335 

A 

3 

0.000790 

-1.12991 

-1.10286 

1.310 

B 

4 

0.000915 

-1.13270 

-1.10369 

1.355 

E 

4 

0.000519 

-1.12511 

-1.10303 

1.335 

A 

4 

0.000790 

-1.12991 

-1.10286 

1.310 

B 

5 

0.000915 

-1.13270 

-1.10369 

1.355 

E 

5 

0.000518 

-1.12509 

-1 . 10303 

1.335 

A 

5 

0.000790 

-1.12991 

-1.10286 

1.310 

B 

true  eigenvalue  =  -1.10263 

*  Values  obtained  by  Goodisman  and  Secrest  (9)  using  their 
A  minimization  scheme. 


■ 

V 

■ 
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calculated  from  the  Temple  formula  and  Weinstein  formula 
respectively.  For  the  Temple  bound  the  value  of  E^ 
was  taken  as  -0.36087  a.u.,  the  exact  eigenvalue  for  the 
first  excited  state  of  £  symmetry  at  2  a.u.  inter- 
nuclear  separation. 

The  second  best  value  of  the  lower  bound  is  obtained 
by  a  modification  of  the  method  of  Fraga  and  Birss,  i.e., 
by  minimizing  a,  as  they  suggest,  but  using  the  Temple 
formula  instead  of  the  Weinstein  formula  to  calculate 
the  lower  bound.  The  lower  bound  obtained  from  the  wave 
function  determined  by  the  ordinary  variational  principle 
is  consistently  the  poorest  value. 

The  convergence  of  the  lower  bound  with  increasing 
number  of  basis  functions  is  quite  rapid.  As  seen  in 
Table  1  there  is  essentially  no  difference  between  the 
results  when  going  from  three  to  five  basis  functions; 
hence  these  values  may  be  considered  to  be  the  best 
lower  bounds  obtainable  with  the  particular  form  of 
trial  function  chosen. 

The  advantage  of  a  self-consistent  iterative  scheme 
such  as  the  one  proposed  here  and  the  scheme  of  Fraga 
and  Birss  is  clearly  shown  when  the  results  obtained 
for  two  basis  functions  are  compared  to  a  calculation 
carried  out  by  Goodisman  and  Secrest  (9) .  The  latter 
calculation  was  based  on  two  terms  of  Eq.  (1)  and  a  two 
dimensional  grid  over  a  and  c^/cq>  to  obtain  the  minimum 
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o£  A.  Their  results  are  included  in  Table  1,  in 
parentheses.  It  is  readily  seen  that  their  values  are 

4 

much  poorer  than  the  true  values  obtained  by  finding  the 
proper  minimum  of  A  by  the  scheme  of  Fraga  and  Birss. 

The  method  of  Goodisman  and  Secrest  also  cannot  be 
extended  easily  to  more  than  two  basis  functions.  Hence 
the  desirability  of  using  a  direct  variational  approach 
(such  as  the  one  of  Fraga  and  Birss,  which  automatically 
determines  the  best  coefficients  and  is  readily  extended 
to  any  number  of  basis  functions)  cannot  be  overemphasized. 

III. 3  Upper  and  Lower  Bounds  to  He  Excited  States 


In  these  calculations,  bas is  functions  of  the  type 
used  by  Davidson  (30)  were  employed.  This  was  done 
because  the  functions  are  relatively  simple  and  also 
the  results  obtained  here  for  the  upper  bounds  could 
be  compared  to  those  obtained  by  Davidson’s  method. 

The  basis  functions  for  the  excited  states  were 


taken  as 


's’n  (n+1)  !  2 


On  ) 


3 

( 2  n  )  ^ 

— C2nsr)n'1eXp(-nsr)Y0j0(6,*) 


n  = 


3 

2 


p’n  (_n+l)  !  2 


— (2n  r)nexpC~n  r)Y-,  M(e,<f>) 
+  n  1 2n+ 1  P  F  P 


n  = 


1.2. 3. 4 

(3) 

2. 3.4.5 

(4) 


where  the  y  are  the  basis  functions  for  the  s  orbital 
As  ,n 

and  the  v  those  for  the  p  orbital.  The  s  and  p 
p  ,n 
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orbitals  were  each  represented  by  four  basis  functions 
and  those  functions  used  for  the  calculations  are  given 
in  Eqs .  (3)  and  (4)  by  the  n  values.  That  is,  the  s 
and  p  orbitals  are  given  by 

4 

d>  =  y  C  V 
s  n=l  s’n  s’n 

5 

d>  =  y  c  v 

p  nt2  p»n  p»n 

Although  the  lx  )  and  (x  }  are  the  same, 

°  s ,  n  p ,  n  ’ 

except  for  the  choice  of  ns  and  rip,  for  all  the  excited 

state  calculations  reported  here,  the  c  and  c  of 

r  *  s ,n  p ,n 

course  differ  for  each  calculation.  They  are  determined 
by  the  direct  application  of  the  method  developed  in 
Chapter  II  for  the  determination  of  the  best  bounds. 
Further  the  n’s  were  varied  such  that  their  optimum 
values  give  the  best  bounds  obtainable  from  a  wave 
function  of  the  type  considered  here.  The  expansion 
coefficients  and  the  optimum  values  of  the  n’s  are  given 
in  Tables  2  and  3  for  the  singlet  and  triplet  states 
respectively.  The  U  or  L  in  parentheses  after  the  state 
designation  in  the  tables  indicate  whether  the  wave 
function  corresponds  to  that  of  an  upper  or  lower  bound 
value  respectively. 

The  values  of  the  upper  and  lower  bounds  calculated 
from  these  wave  functions  are  given  in  Table  4  along 
with  the  experimental  energies  and  the  values  obtained 


Table  2 


Singlet  Excited  State  Wave  Functions  -  Bound  Method 


n 

n 

c 

n 

n 

n 

cR 

1P 

Is  2p 

CL) 

1 

2.19 

0 .87274 

2 

0.368 

-0.02899 

2 

2.19 

0 .13117 

3 

0 . 368 

0.02990 

3 

2.19 

0.01544 

4 

0.368 

0.99266 

4 

2.19 

0.00195 

5 

0 . 368 

0.00057 

1P 

Is  3p 

(U) 

1 

2.19 

0 .87274 

2 

0.364 

-0.04185 

2 

2.19 

0 . 13117 

3 

0.364 

0.00600 

3 

2.19 

0.01544 

4 

0.364 

1.03134 

4 

2.19 

0.00195 

5 

0.364 

-0.00699 

1P 

Is  3p 

(L) 

1 

2.13 

0.90986 

2 

0.269 

0.01350 

2 

2.13 

0.09619 

3 

0.269 

-0.02806 

3 

2.13 

0.00800 

4 

0.269 

-0.94830 

4 

2.13 

0.00068 

5 

0.269 

1.86615 

1P 

ls4p 

(L) 

1 

2 .13 

0.90986 

2 

0 . 180 

-0.07002 

2 

2.13 

0.09619 

3 

0.180 

1.10376 

3 

2.13 

0.00800 

4 

0.180 

-3.40898 

4 

2.13 

0.00068 

5 

0 .180 

3.30019 

• 

• 

Table  3 


Triplet  Excited  State 

Wave 

Functions  • 

-  Bound  Method 

<\> 

s 

<J> 

P 

n 

n 

c 

n 

n 

n 

cR 

3p 

Is  2p 

(L) 

1 

1.76 

1.21140 

2 

0.377 

-0.07442 

2 

1.76 

-0 .28595 

3 

0.377 

0.23490 

3 

1.76 

0.05236 

4 

0.377 

0 . 84289 

4 

1.76 

-0.00705 

5 

0.377 

-0.00950 

3p 

ls3p 

(U) 

1 

1.76 

1.21140 

2 

0.376 

-0.10792 

2 

1.76 

-0.28595 

3 

0.376 

0 . 26166 

3 

1.76 

0.05237 

4 

0.376 

0.78938 

4 

1.76 

-0.00706 

5 

0 .376 

0.04704 

3P 

Is  3p 

(L) 

1 

1.85 

1.12406 

2 

0 .276 

0.03428 

2 

1.85 

-0.15784 

3 

0.276 

-0.17941 

3 

1.85 

0.01740 

4 

0.276 

-0.60847 

4 

1.85 

-0.00153 

5 

0.276 

0.16616 

3p 

ls4p 

CL) 

1 

1.86 

1.11500 

2 

0.186 

-0.12947 

2 

1.86 

-0.14535 

3 

0.186 

1.27357 

3 

1.86 

0.01489 

4 

0 . 186 

-3.73863 

4 

1.86 

-0.00122 

5 

0.186 

3.45304 
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by  using  Davidson’s  method  with  the  basis  set  discussed 

above  (E  ...  -. )  ,  and  with  a  larger  basis  set 

(E ’  .  , .  ,).  A  discussion  of  Davidson’s  method  is  in 

v  variational 

order  now  before  continuing. 


Outline  of  Davidson’s  Method 

For  the  two  electron  problem  considered  here  the 
trial  function  may  be  written  as 


1.3  = 


*■’'  =  1/^  [*S(D*  (2)  *  *  (l)*s(2)]  (7) 


where  <J>  and  4>  are  as  given  in  Eqs .  (5)  and  (6)  ; 
s  p 

<p  =  c  x  and  4>  =  c  x  •  Following  the  usual  notation 

s  s  s  p  p  p 


H  =  H^l)  +  H1  (2)  +  l/r12 


(8) 


H-  .  =  /  x  •  x^  d- 
~ij  J  i  1  J 


(9) 


S.  .  =  /  txt  dx 

~ij  J  iAj 


(10) 


Js  =■  /  X  Cl)X5(2)csrjlFxs(2)XpCUdT1dT2 


-!■?+: 


(ID 


Ks  =■  /  X  Cl)xsCl)csr-12^xs(2)Xp(2)dT1dT 


(12) 

with  similar  definitions  for  J  and  K  .  Minimizing  the 

~P  _  p  a 

average  energy,  <$|H|<f>>,  with  respect  to  the  coefficients 
yields  Roothaan’s  (33)  pseudoeigenvalue  equations 


■ 


. 
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(H '  +  e  S  )c  =  E  (§  ±  §  c  c+S  )c 

^..S  p„SS  S  _;SS  Isp  P  p^ps  s 


(13) 


(H '  +  e  S  )c  -  E  ±  S  c  c  S  )c 

,p  s.pp'  p  _pp  .ps  s  s.sp;  p 


C14) 


—V 

where  c  and  c  are  such  that  c  -S.  .c  •  =  1,  and 
s  p  1.111  * 


H  ’ 
-P 


H  +  J  ±  (K  +  H  c  cfS  +  S  c  c+H  ) 
~pp  ~s  .s  .ps  s  s.sp  „ps  s  S~Sp' 


and 


-v  f  -> 

e  =  c  -  H .  .c  .  . 

S  1.111 


(15) 


(16) 


Since  xg  and  Xp  are  of  different  symmetry,  Ssp  =  0 
and  Eqs .  (13)  and  (14)  simplify  to 


H'c  =  A  S  c 
.s  s  s.ss  s 


(17) 


also 


and 


H’c  =  A  S  c 
~P  P  P~PP  P 


E  =  A  +  e  =  A  +  A  -  c  +  (J  ±  K  )c 
s  p  s  p  p.s  ~s  p 


(18) 


(19) 


H f  =  H  +  J  1  K 
.p  .pp  ~s  ~s 


(20) 


Now  suppose  the  eigenvalues  of  Eq.  (17)  are 


A  ,  _<  A  ?  <  .  .  .  A  M  with  the  corresponding  eigenvectors 

S  j  1.  S  j  Z  S  ^  IN 

-y  ~y  *y 

c  ,,  c  ?  ...  c  N  and  similarly  for  Eq.  (18). 

S  ,  1  S  ,  Z  S  y  IN 

Davidson  points  out  that  any  of  the  c  . ,  say  c  can 

S  j  1  S  j  Ct 

“V 

be  used  to  form  HT  and  any  of  the  c  . ,  say  c  „  can  be 

-P  7  p ,  i  *  7  p ,  8 

used  to  form  HT  and  hence  at  self-consistency  Eqs.  (17) 

~  s 


■ 


V 

. 

. 
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and  (18)  will  be  self-consistent  with  respect  to  c  (a, 3) 

S  -y  CL 


and  c  fa  J)  .  Davidson  proyes  that 


<t  HU.>  =  6.  . 
k  1  kj 


(21) 


where 


and 


K  =  v(2)  1  v(l)^(2)] 


p  ,k 


p,k 


<f>s  =  CIC1,B)XS.  ^p  k  =  cp  k^’^^p* 


Hence  the  configurations  are  orthogonal  and  non-inter¬ 
acting.  By  using  the  method  of  MacDonald  (34)  which  gives 
upper  bounds  to  eigenvalues  from  the  solution  of  the 
secular  equation  when  the  wave  function  is  expanded 
in  a  truncated  set  of  complete  functions,  one  may  write 


‘sCWn)  +  xp,b(1’6)  i  eb 
th 


(22) 


where  is  the  true  g  eigenvalue  of  H.  Thus  one  is 
provided  with  a  variational  principle  for  excited  states 
in  this  simple  case  which  insures  an  upper  bound  to  the 
true  eigenvalue.  This  method  provides  a  convenient 
comparison  between  upper  bounds  calculated  by  an  ordinary 
variational  treatment  and  the  bounds  calculated  by  the 
method  proposed  in  this  thesis. 

The  method  just  outlined  provides  the  means  of 
obtaining  the  values  listed  in  Table  4  in  the  fifth 


and  sixth  columns.  The  values  in  the  fifth  column  were 


Nv 


- 


. 
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obtained  by  direct  application  of  Davidson's  method  . 
using  the  basis  sets  described  above  and  optimizing  the 
n ’ s .  The  wave  functions  corresponding  to  these  energies 
are  given  in  Tables  5  and  6.  It  should  be  noted  that 
some  of  the  vectors  given  in  Tables  5  and  6  as  well 
as  some  of  those  in  Tables  2  and  3  have  coefficients 
which  are  greater  than  unity.  These  vectors  are, 
however,  normalized,  as  they  must  be;  the  reason  for 
the  strange  coefficients  is  the  form  of  the  functions 
chosen,  i.e.,  all  the  basis  functions  of  the  same 
symmetry  have  a  common  n .  Usually  each  basis  function 
has  its  own  orbital  exponent  which  is  optimized  and  the 
coefficients  then  are  all  close  to  or  less  than  unity. 

The  sixth  column  of  Table  4  contains  the  results 
obtained  by  Davidson  using  larger  basis  sets;  these 
results  should  be  close  to  the  Hartree-Fock  limit. 

It  should  be  noted  that  the  bounds  calculated 

give  quite  good  limits  on  the  true  eigenvalues.  The 

difference  between  the  upper  and  lower  bound  values 

being,  in  the  worst  cases,  only  .005  a.u.  (~'3  kcal.). 

For  the  lowest  states  of  a  given  symmetry  the  Upper 

Bound  value  and  the  E  .  . .  value  are  identical,  since 

variational  * 

the  upper  bound  equation  reduces  to  the  ordinary  varia¬ 
tional  principle  in  these  cases . 

The  fact  that  the  upper  bounds  given  by  the 

Davidson  method  fE  ...  of  Table  4)  are  better 

v  variational 


i 

- 

' 

. 

10® 


. 

*  •  ••  . .  r  01  ,a*  I' 


1 

2 

3 

4 

1 

2 

3 

4 

1 

2 

3 

4 


Table  5 


State  Upper  Bound  Wave  Functions  -  Davidson' 


<f>s  •  <t»p 


T) 

c 

n 

n 

c 

n 

n 

XP 

ls2p 

(A) 

1.64 

1.34383 

2 

0.624 

0.46829 

1.64 

-0.50259 

3 

0.624 

0.48996 

1.64 

0.14015 

4 

0.624 

-0.03214 

1.64 

-0.02830 

5 

0.624 

0.12795 

1P 

Is  3p 

(A) 

1 . 80 

1.17168 

2 

0.355 

-1.25864 

1.80 

-0.22546 

3 

0 . 355 

1.13821 

1 . 80 

0.03362 

4 

0.355 

0  .  76193 

1.80 

-0.00449 

5 

0.355 

-0.03427 

1P 

ls4p 

(A) 

1.88 

1.10042 

2 

0.257 

1.97148 

1.88 

-0.12536 

3 

0.257 

-4.84864 

1.88 

0.01121 

4 

0.257 

2.61793 

1.88 

-0.00113 

5 

0.257 

0.96264 

Method 


n 

1 

2 

3 

4 

1 

2 

3 

4 

1 

2 

3 


Table  6 


State  Upper  Bound  Wave  Functions  -  Davidson’s  Method 


^  s  ‘f’p 


n 

cn 

n 

n 

cn 

• 

3P 

Is  2p 

(A) 

1 . 76 

1.20304 

2 

0.710 

0.68640 

1.76 

-0.27858 

3 

0.710 

-0.10798 

1.76 

0.05158 

4 

0.710 

0 . 50040 

1.76 

0.00094 

5 

0.710 

0.00028 

3P 

ls3p 

(A) 

1 . 79 

1.17834 

2 

0.437 

-1.17537 

1.79 

-0.23779 

3 

0.437 

1.42767 

1.79 

0.03893 

4 

0.437 

-0.79388 

1.79 

-0.00328 

5 

0.437 

1.17645 

3P 

ls4p 

(A) 

1.85 

1.12427 

2 

0.257 

2.28960 

1.85 

-0.15873 

3 

0.257 

-6.40743 

1.85 

0.01798 

4 

0 .257 

5.01656 

1.85 

-0.00101 

5 

0.257 

-0.21458 

■ 
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than  those  of  the  method  proposed  in  this  thesis  for 
1  3 

the  P[ls3p)  and  PCls3p)  states  is  explained  by  a 
simple  argument.  In  the  special  case  treated  by  Dayidson, 
his  method  is  equivalent  to  a  variational  principle 
which  determines  the  function  which  gives  a  stationary 
value  to  E,  i.e.,  it  satisfies  the  equation,  6E  =  0. 

This  is  the  same  form  as  the  ground  state  variational 
principle  and  it  is  only  because  of  the  nature  of  the 
trial  functions  chosen  and  the  form  of  the  problem  which 
allows  this  equation  to  be  valid  for  the  determination  of 
excited  state  functions  also.  As  was  shown  previously 
the  method  proposed  here  gives  an  equation  of  the  form 
<5E  +  w  6  A  =  0. 

The  latter  equation  may  be  considered  as  a  constrained 
variational  problem  with  the  former  being  the  unconstrained 
problem.  These  problems  are  well  known  and  it  is  a  fact 
that  the  constrained  problem  cannot  have  a  minimum  which 
is  as  low  as  the  minimum  for  the  unconstrained  problem. 
Hence  Davidson's  values  are  better  upper  bounds.  How¬ 
ever  in  the  general  case,  when  Davidson's  method  is 
not  applicable,  the  other  method  is  still  valid  and 
will  yield  values  for  an  upper  bound  which  would  be 
otherwise  unavailable. 

I I I . 4  He  Excited  State  Wave  Functions  by  e2/A  Method 


In  section  1.3  a  method  was  proposed  for  the 


I 


V 


S. 


. 
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determination  of  an  approximate  wave  function  for  excited 
states  which  has  the  maximum  overlap  with  the  true 
eigenfunction.  The  method  is  based  on  the  minimization 
of  £ 2 / A  as  described.  In  Chapter  II  the  method  was 
developed  in  terms  of  orbital  theory  and  it  is  that 
formalism  which  can  be  directly  applied  here  to  yield 
the  excited  state  wave  function.  The  calculations  were 
carried  out  for  four  states  of  He:  '*'P(ls2p),  ‘*'P(ls3p), 
^P(ls2p)  and  ^P(ls3p)  states.  The  basis  set  used  is  the 
same  as  that  used  for  the  calculations  described  in  the 
last  section. 

For  these  calculations  one  should  know  the  in 

order  to  calculate  e.  In  all  cases  considered  here  the 

were  approximated  by  using  the  experimental  energies 

of  the  appropriate  states.  It  should  be  pointed  out 

that  if  one  did  not  know  the  W,  or  have  any  experimental 

k 

information  that  a  double  minimization  procedure  could 
be  devised  to  obtain  the  best  approximation  to  W^. 
automatically.  That  is,  one  could  find  the  minimum  of 
.£ 2 / A  for  a  series  of  values  of  then  considering  the 

minimum  values  as  a  function  of  one  would  choose  the 

smallest  of  these  minimum  values  and  thus  obtain  the 
wave  function  and  an  approximation  to  the  corresponding 


V 


In  Table  7  the  wave  functions  determined  by  the 


i  1 
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Table  7 


State  Wave  Functions 


e2/A 


Minimization  Method 


n 

cn 

n 

n 

cn 

1P 

Is  2p 

(B) 

2 .25 

0 .83850 

2 

0.668 

0.41858 

2.25 

0 . 16110 

3 

0.668 

0.44990 

2.25 

0.02780 

4 

0.668 

0.02159 

2.25 

-0.00146 

5 

0.668 

0 .17857 

3 

Is  2p 

(B) 

1.86 

1.11008 

2 

0 . 700 

0.71746 

1.86 

-0.14467 

3 

0 . 700 

-0.10578 

1.86 

0.01676 

4 

0.700 

0.46560 

1.86 

0.00726 

5 

0.700 

-0.00067 

1P 

Is  3p 

(B) 

1.68 

1.29512 

2 

0 . 354 

-1.31116 

1.68 

-0.42127 

3 

0.354 

1.29494 

1.68 

0 . 10502 

4 

0.354 

0.62949 

1.68 

-0.01848 

5 

0.354 

-0.00507 

3P 

ls3p 

(B) 

2.45 

0.73860 

2 

0.436 

-1.20234 

2.45 

0.23844 

3 

0.436 

1.48308 

2.45 

0.03785 

4 

0.436 

-0.78487 

2.45 

0.03638 

5 

0.436 

1.14108 

. 

■ 
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£ 2 / A  minimization  meth 
presented  in  Chapter  I 
the  functions  obtained 
overall  wave  function, 
section  by  calculating 
operators  important  in 


od  are  presented.  The  arguments 
for  this  method  maintained  that 
by  this  method  should  be  the  best 
This  is  investigated  in  the  next 
some  expectation  values  of 
physical  properties. 


III. 5  Calculation  of  Some  Expectation  Values 


In  order  to  demonstrate  the  validity  of  the  thesis 
that  the  wave  functions  obtained  by  the  minimization 
of  e2/A  should  be  the  best  functions  over  the  whole  of 
configuration  space,  some  expectation  values  have  been 
calculated  for  these  functions. 

If  the  trial  function  is  designated  as  $,  the 
expectation  values  calculated  are  given  as 


E  =  <  $  I  H  I  $> 


<l/r>  =  l  <$|  l  1/r. 

z  i=l  1 


$> 


<r>  = 


,  2 

2  <*I.I 

z  i=l 


ril*> 


(23) 

(24) 

(25) 


<r2>  =  \  <$i  y 
2  i=l 


r ?  I  $> 


(26) 


.  th 


where  r^  is  the  distance  of  the  i  electron  from  the 
nucleus.  In  Table  8  a  comparison  is  made  of  the 


expectation  values  calculated  from  the  P(ls2p)  and 


$  i 
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3 

P(ls3p)  wave  functions  of  Table  7  with  the  expectation 
values  calculated  from  variational  wave  functions 
obtained  by  other  methods.  The  other  wave  functions  are 
those  designated  as: 

1)  -  the  wave  functions  of  Tables  2  and  3  which  are 
there  designated  by  (L)  and  described  in  detail  in 
section  1 1 1 . 3 , 

2)  E  -  the  appropriate  wave  functions  of  Tables  5  and  6 
which  are,  in  these  cases  of  the  lowest  state  of  a 
given  symmetry,  equivalent  to  the  wave  functions 
determined  by  the  ordinary  variational  principle. 

3)  "Exact”  -  the  wave  functions  determined  by  Schiff 
et .  al .  (35)  using  560  terms  of  a  Pekeris  type 
expansion.  The  latter  are  extremely  accurate  functions 
and  for  all  practical  purposes  may  be  considered 

to  be  exact  -  hence  the  designation. 

Comparing  the  values  in  the  first  three  columns 
of  Table  8  with  those  of  the  fourth  column,  it  is  readily 
seen  that  the  e2/A  wave  functions  give  the  best  overall 
accuracy;  indeed,  the  only  value  predicted  better  by 
the  variational  principle  is  that  of  the  energy  and  the 
difference  there  is  in  the  fifth  decimal  place. 

A  similar  comparison  is  given  for  the  ^PCls3p)  and 
^P(ls3p)  states  of  He  in  Table  9.  For  these  states 
however  the  ordinary  variational  principle  cannot  be 
used  and  hence  the  method  of  Davidson  described  in  the 
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Table  8 


Comparison 

of  Expectation  Values 

Calculated 

from 

-  _  \ 

Various 

Wave  Functions:  ls2p 

1?  Is  2p 

bl 

E 

e  2/ A 

"Exact'’ 

-E 

2.07204 

2.12246 

2.12243 

2.12384 

<l/r> 

1.0457 

1.1194 

1.1229 

1.12318 

<r> 

6.134 

2.9469 

2.9223 

2.91068 

2 

<r  > 

83.87 

16.175 

15.878 

15.76565 

3P  ls2p 

bl 

E 

£  2  /  A 

"Exact" 

-E 

2.07537 

2.13134 

2.13132 

2.13316 

<l/r> 

1.0490 

1.1309 

1.1335 

1.13324 

<r> 

6.185 

2.7004 

2.6813 

2.67396 

<r2> 

75.54 

13.417 

13.245 

13.21174 

Table  9 


 Comparison 

of  Expectation  Values 

Calculated 

from 

Various 

Wave  Functions  :  ls3p 

1P 

ls3p 

bl 

BU 

E 

e  2/A 

’'Exact” 

-E 

2.03844 

2.07062 

2.05473 

2.05472 

2.05515 

<l/r> 

1.0448 

1.0229 

1.0547 

1.0528 

1.05497 

<r> 

6.122 

12.055 

6.720 

6.695 

6.67954 

2 

<r  > 

86.77 

296.20 

92.81 

92.16 

91.8727 

3p 

ls3p 

bl 

BU 

E 

e2/A 

"Exact” 

-E 

2.03999 

2.07374 

2.05750 

2.05745 

2.05808 

A 

M 

\ 

H 

V 

1.0474 

1.0239 

1.0570 

1.0597 

1.05802 

<r> 

6.302 

11.574 

6.364 

6.332 

6.32113 

<r2> 

88.48 

268.22 

82.96 

82.15 

82.1098 
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last  section  applies  to  the  determination  of  the  wave 
functions  in  Tables  5  and  6.  It  is  the  appropriate  wave 
functions  of  these  tables  that  were  used  to  calculate 
the  expectation  values  given  under  the  designation  E 
in  Table  9.  The  other  designations  of  Table  9  have  the 
same  significance  as  in  Table  8  with  the  addition  that 
denotes  the  results  obtained  from  the  appropriate 
wave  functions  given  in  Tables  2  and  3.  Again  the  method 
involving  the  minimization  of  e2/A  gives  the  best  overall 
agreement  with  the  "Exact”  values  of  the  operators 
considered. 

Tables  8  and  9  demonstrate  a  rather  important  result: 
an  approximate  wave  function,  is  quite  sensitive  to  the 

criterion  which  is  used  to  determine  it,  and  although 
it  may  predict  the  property  very  well  which  was  used 
to  judge  its  "goodness",  other  properties  can  be  very 
inaccurate.  This  is  particularly  the  case  of  the  func¬ 
tions  which  give  the  upper  and  lower  bound  values. 

Although  these  functions  give  very  good  bounds,  as  shown 
in  Table  4,  they  are  practically  useless  for  determining 
any  other  property.  To  a  lesser  degree  the  same  criticism 
of  the  variational  principle  can  be  made;  in  trying  to 
obtain  the  lowest  possible  value  of  E,  a  rather 
significant  sacrifice  is  made  in  the  wave  function  with 
respect  to  its  ability  to  give  reliable  values  for 
other  properties. 
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The  minimization  of  ,e2/a  gives  a  general  yariational 
principle,  which  is  applicable  to  the  ground  or  excited 
states  regardless  of  symmetry,  for  the  determination 
of  the  best  trial  wave  function  <j> . 

III. 6  Calculations  of  Lower  Bound  to  He  Ground  State 


As  was  mentioned  at  the  end  of  Chapter  I,  there 
exist  a  few  non- variational  approaches  to  determining 
lower  bounds  of  eigenvalues.  Of  these  methods  the  one 
of  Lowdin,  which  is  based  on  perturbation  theory,  has 
received  a  great  deal  of  attention.  It  is  the  purpose 
of  this  section  to  make  a  direct  comparison  between  the 
method  proposed  in  this  thesis  and  the  method  of  Lowdin, 
by  considering  the  results  of  the  two  methods  for  the 
case  of  the  ground  state  of  the  He  atom. 

The  trial  functions  used  for  the  method  proposed 
here  are  of  the  general  type 


where  the  C 
variational 
the  maximum 


*  =  I  c„_x 


£mn 


£mnA£mn 


£mn 


are  coefficients  to  be 


calculation  and  will  be  s 


value  of  the  lower  bound. 


(27) 

determined  by  the 
uch  that  $  gives 


•£mn 


XT  £.m  n  -s/2 
=  N„  s  t  u  e 
£mn 


C28) 


where  N 

1  =  ko2 


mn 


is  a  normalization  factor,  s  =  k(r^  +  r?), 
r-^)  and  u  =  kr-^*  m  an  even  integer  and 
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k  is  a  scaling  factor.  The  form  of  the  Hamiltonian  and 
all  the  necessary  integrals  are  discussed  in  Appendix 
IV.  For  the  calculations  presented  here  the  values  of 
Hmn  taken  in  the  summation  for  the  three  trial  functions 
considered  are  presented  in  Table  10.  The  latter  are 
also  the  same  functions  used  for  the  "inner  projection" 
space  in  the  Lowdin  method  calculations  presented  here. 

A  derivation  of  Lowdin' s  method  is  given  in  Appendix 
III  and  the  necessary  integrals  are  given  in  Appendix  IV. 

In  Table  11  the  results  of  the  two  methods  are  com¬ 
pared  with  the  results  of  the  Temple  method  for  three 
basis  sets.  The  basis  functions  listed  in  column  one 
of  Table  10  were  used  for  the  6-function  basis  set. 

The  functions  of  column  two  were  added  for  the  10-function 
calculations  and  all  of  the  functions  shown  to  Table  10 
were  used  for  the  18-function  calculations.  In  Table  11 

* 

the  numbers  in  parentheses  are  the  optimum  values  of 
the  scaling  factor  k. 

The  results  show  that  both  the  present  method  and 
Lowdin' s  method  give  a  considerably  better  lower  bound 
than  does  the  usual  Temple  method.  The  difference 
between  the  present  method  and  that  of  Lowdin  for  18 
basis  functions  is  not  very  great  and  from  the  trend 
of  Table  11  larger  basis  sets  would  probably  give 
almost  identical  results.  This  is  a  rather  important 
point  since  the  present  method  can  be  extended  to 
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Table  10 


Basis  Functions  for  Calculations 


m  n 
0  0  0 
0  0  1 
0  2  0 
10  0 
2  0  0 
0  0  2 


m  n 
10  1 
0  2  1 
0  0  3 
0  2  2 


m  n 
12  0 

3  0  0 

0  2  4 

0  0  4 

0  0  5 

0  2  3 

2  2  0 

4  0  0 


Table  11 


Comparison  of  Lower  Bound  Values 


Basis  Functions  Present  Method 


t 


Lowdin  Method*  Temple  Method 


t 


6 

10 

18 


-2.9147  (3.34) 

-2.9080  (3.91) 

-2.9041  (4.02) 


-2.9098  (3.46) 

-2.9060  (3.50) 

-2.9039  (3.92) 


-2.9260  (3.64) 

-2.9132  (3.54) 

-2.9055  (3.85) 


t  E^  =  -2.1460  was  used  as  the  lower  bound  to  the 

true  eigenvalue  of  the  2  state 
of  He  (37) . 

*  £  =  -2.9037  was  used  as  the  upper  bound  to  the 

true  eigenvalue  of  the  1  state 
of  He  (21) . 
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handle  many  electron  atoms  without  much  difficulty, 
whereas  Lowdin's  method  has  the  simple  form  Csee 
Appendix  III)  discussed  here  only  for  two  electron 
systems,  since  the  perturbation  term  is  simply  V  =  1/r-^* 
For  many  electron  atoms  the  perturbation  is  much  more 
complicated  and  the  integrals  which  would  arise  as  a 
result  of  this  are  quite  intractable  at  present. 

III. 7  Summary  and  Conclusions 

It  has  been  shown  here  that  it  is  possible  to 
construct  variational  methods  for  the  determination  of 
upper  and  lower  bounds,  and  for  the  determination  of  the 
best  approximate  wave  function  in  the  sense  of  the  maximum 
overlap  with  the  true  eigenfunction.  The  results  have 
shown  that  these  methods  are  not  only  workable  but  that 
they  yield  quite  good  values. 

For  the  bounds  determined  here,  the  associated  wave 
functions  have  proved  to  be  useless  for  the  calculation 
of  physical  properties.  In  the  Lowdin  method  for  lower 
bounds  the  corresponding  wave  function  cannot  be  written 
in  a  tractable  form.  However  it  would  most  certainly 
suffer  from  the  same  characteristic  lack  of  utility,  which 
seems  to  be  the  fate  of  wave  functions  associated  with  a 
bound  value.  Hence  the  variation  method  developed  here 
for  bounds  is  quite  useful  in  determining  bounds  but 
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should  not  in  any  way  be  considered  as  a  way  to  determine 
a  valid  approximation  to  an  eigenfunction.  The  reason 
is  that  the  bound  formulations  add  contributions  from 
other  configurations  to  the  trial  wave  functions  which  do 
not  help  in  the  description  of  the  wave  function  associated 
with  the  state  under  consideration.  This  is  most  clearly 
seen  in  Lowdin's  formulation  for  lower  bounds,  where 
contributions  to  the  reference  function  are  added  from 
functions  projected  from  another  subspace.  It  is  also 
true  for  the  method  proposed  here.  The  validity  of  this 
statement  can  be  attested  to,  numerically,  by  the  wave 
functions  presented  here  for  the  bounds.  To  obtain  good 
bounds  one  must  add  in  parts  of  function  space  which 
are  not  suited  to  the  description  of  the  eigenfunction. 

The  ordinary  variational  principle  is  a  special  case  of 
the  bound  method  developed  here.  However  it  is  a  special 
case  which  fortunately  does  not  suffer  from  the  draw¬ 
back  just  discussed.  This  is  due  to  the  fact  that  the 
function  satisfies  the  equation  <5E  =  0,  rather  than 
6E  +  co6 A  =  0.  It  is  the  latter  which  is  responsible  for 
adding  in  the  terms  which  are  undesirable  from  the 
point  of  view  of  obtaining  a  good  wave  function;  in 
particular,  it  is  the  quantity  w  which  is  responsible. 

This  can  be  seen  by  comparing  the  wave  functions 
obtained  by  the  bound  method  with  those  obtained  by 
the  e2/A  method.  In  the  former  the  oo  is  of  such  a  form 
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that  it  adds  in  undesirable  contributions,  but  in  the 
latter  the  w  is  of  the  correct  form  to  obtain  the  best 
function.  Although  the  ordinary  variational  method  is 
thus  a  reasonable  method  for  obtaining  approximate  wave 
functions,  this  is  only  of  a  secondary  nature.  The 
primary  utility  of  the  variational  principle  for  lowest 
states  should  be  recognized  as  a  method  for  determining 
the  best  upper  bound  to  the  lowest  state  of  a  given 
symmetry . 

For  the  determination  of  the  best  approximate  wave 
function  a  variational  method  based  on  the  minimization 
of  e2/A  has  been  developed.  The  resultant  wave  function 
has  the  maximum  overlap  with  the  true  eigenfunction . 
Unlike  the  ordinary  variational  principle  (which  really 
gives  the  best  upper  bound  for  the  ground  state  and  not 
the  best  approximation  to  the  eigenfunction) ,  the  method 
proposed  here  gives  the  wave  function  which  has  the 
least  mean  square  deviation  from  the  true  eigenfunction 
over  the  whole  of  configuration  space.  The  results  on 
He  indeed  show  that  this  method  is  superior  to  the 
ordinary  variational  principle  for  the  determination  of 
wave  functions.  It  is  furthermore  general  and  can  be 
used  to  determine  an  approximate  wave  function  for  any 
state . 

In  quantum  chemistry  one  is  not  only  interested  in 
the  energy  of  a  system  but  also  in  calculating  physical 
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and  chemical  properties.  Since  the  method  proposed  here 
yields  wave  functions  which  are  capable  of  giving  values 
of  physical  properties  which  are  superior  to  those  of 
the  ordinary  variational  methodjit  is  hoped  that  this 
may  constitute  a  contribution  to  the  better  prediction 
of  phenomena  of  chemical  interest. 
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Appendix  I 


According  to  Eq .  (2.2)  E  may  be  written  as 
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Now  considering  each  term  separately  one  has 
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The  evaluation  of  the  matrix  components  of  the  Hamiltonian 

with  respect  to  Slater  determinants,  as  given  in  the 

equations  above,  can  be  found  in  standard  texts  on  quantum 

mechanics  (23),  however  the  notation  here  is  quite  different. 

The  l  is  to  be  interpreted  in  the  following  way.  The 
i  Aa 

function  $  consists  of  a  set  of  orbitals 
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This  set  of  orbitals  may  be  decomposed  into  n  subsets 
according  to  the  symmetry  species  and  subspecies  of  the 
orbitals,  hence  one  has  the  group  of  subsets 
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where  d(Aa)  is  the  dimension  of  the  subspace  spanned  by 
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the  subset  elements  d> .  .  The  summation  is  then  to  be 
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carried  out  over  all  subsets  and  over  all  orbitals  in  a 
subset.  This  is  obviously  equivalent  to  a  summation  over 
all  orbitals  in  the  original  set. 

The  following  notation  will  also  be  introduced: 

£  =  l .  Thus 

iXa  t 


y  fAHXa  e  y  f?H^. 

.  f-  11  ft  1  1 

l  X  a 


i 


Note  that  the  fX  are  not  summed  over  the  index  a  and  hence 

i 

the  shorthand  notation  for  f^  cannot  be  used.  In  all  such 
cases  where  the  complete  set  of  indices  i_s  not  summed  over 
the  subset  of  indices  which  i_s_  summed  over  will  be  designated 
by  explicitly  affixing  these  indices  to  the  quantities 
involved . 

Using  this  new  notation, 
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i 


imm 


lm 


where  the  and  are  as  follows 

li  limm 


H  H  d)  -*■> 

lm  l 1  1 rm 


(1.4) 


(1.5) 


lmnp 


2 


j  N 

lmnp ; 1  t ; lmnp 


imnp ; 2 


N 

4-  ;  imnp 


(1.6) 


where 


N  —  <  (i  I  N  A 

i ; imnp  Yi 1  t ;np 1 Ym 


If 


<f>J(P)  <J>^(  p)  [ !/ r p T ]  4^(t)  <J>  J(T)  dVpdVT 


ri.7) 


■ 
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and  N.  is  the  exchange  counterpart  of  (1.7),  i.  e., 

1 ; lmnp  &  1  v  * 

(1.7)  with  the  electron  labels  of  the  second  and  fourth 


orbitals  interchanged.  The  nj 


and  n „  are  constants 
i mnp ; 2 


lmnp ;1 

which  depend  upon  the  occupancies  of  the  various  orbitals 
involved  and  upon  the  multiplicity  and  orbital  angular 
momentum  of  the  multiplet  under  consideration. 
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Appendix  II 


From  Eq .  (2.5) 


H2  =  l  +  2  l  l  +  l  l  J_  J_  (II.1) 

TTO  TT  p  <  T  pi  TT<a  p<T  7T  O  pT 


Now  splitting  the  first  and  third  terms  into  two  new  terms 


each , 


9  9  H 

h  =  y  hz  +  y  hh  +2  y  y  — 

“  TT  H  TT  O  L  L  y 

IT  Of  TT  7TP<XPT 


l  l 


1  1 


XT  <0  r 


TT  0 


^  ^  r  r 
XT<G  p  <  T  TT  O  pi 


(II. 2) 


Now  splitting  up  the  third  and  fifth  summations  of  (II. 2) 
one  obtains  finally, 


h2  =  y  h2  +  y  hh  +  2  y  - 

L  xr  5  xt  a  u  r 

XT  O  jr  XT  XT  <  T  TT  T 


H  H 

'  +  2  I  I  ^ 


TT  p  <T  P  T 

P  7^  TT 


i  1 


i  4-*  z  i  4^  +  i  i  r  r 

XT  <  O'  r  TT  <  O  0<T  XT  O  OT  TT  <  O  P<X  XT  o  px 

TT^X  TT^pO^X  71 

a^p 


r2  . 


This  arrangement  of  H  is  more  convenient  for  evaluating 

2  2 
<$|H  |  3>>  and  obtaining  unique  terms.  Each  term  o£  H 

will  now  be  evaluated  with  respect  to  $. 


y  < $  | h^  |  $>  =  y  f.xi-i-  =  y  f .y^ 

L  1  XT  1  ±  1A11  X  1A1 


XT 


i 


i 


(II. 4) 


' 

- 
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tt 


Tf  a 


=  -2  l 

f  XfyH->->H 

l  m  im  : 

im 

i^m 

H 

y  <$|  71 

L  1  r 

1  $>  -  1 

1  -> 

IT  <  T  TT  X 


I  fXfyH^->H+->  (1 1. 5) 
1  m  1 1  mm  v  J 
lm 


lm 


£vf 

l  m  lm 


(II. 6) 


l  <$ 

tt  <a 


2  I  l  <* 

TT  p  <T 
P^TT 


4-1  $>  =  l  fXfyQ 
2  1  -*->  i  mx 


tt  a 


H 


TT 


P  T 


im 


lm 


(II. 7) 


$>  =  -  I  f*fvfu 

i  m  n  im  minn 

lmn 

ifm 


+  3  T  fXfyfvH^N-^->->  (II.  8) 
^  i  m  n  n  mmnn  v  J 

lmn 


l  y  <$|  — — — |  $> 

L  L  I  ^  IT 

tt  <a  a  <x  ttcj  ax 

TT^X 


7  fXfyfv(l>++  +  Y++->) 
l  m  nv  lmn  imrr 

imn 


(II. 9) 


v  v  ,i  1  1 

l  l  <$l^ £ 

TT<ap<X  TTG  PX 


TT^P  a^x 
a^p 


$>  =  -  V  fXfyfvfC 

l  m  n  p  nmn  nmpp 

imnp 

m^n 


++++  3  imnp 1  pmni  nmpiJ 
imnp 


l  (II  .10) 

iimm  nnpp  v  J 

imnp 


where  and  Hi-*  are  as  defined  in  Appendix  I  and  where 

imnp  im 


Pi"^  =  2pi"^  4  P  4  -t-*-  -  p ~P ,  -*-»■ 
im  rim;i  +;im  rim;2  +;im 


(11.41) 


=  2q^+  j  Q  „Q ,  -} 

xim  nim;lxt;im  nim;2x4-;i 


im 


(II  .12) 


imn  imn;l  f;imn  imn ; 2  + ; imn 


(II  .13) 


' 


. 


V.  <t 
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lmn 


2y->+->  i  Y 
' lmn ; I  t ; 


imn  '  lmn  ;  2  4- ;  lmn 


(II  .14) 


The  explicit  forms  of  the  integrals  involved  in  the  above 
expressions  are  given  below. 


Xim  =  I  [H%;(p)]  [H^(p)]dV 
1  ^ 


(II  .d5) 


P  =  //  4>*  CpD  4>*^Cp3  C - )  4>*  (t)  [ H 4>->  C T )  ]  dV  dV 

+  ;im  JJ  K  J  l  vr  '  J  L  p  T 

’  l  p  t  .  m 


1 1 

■t ;  lm 


d  ,2 


=  //  4>J(p)  <f>j(p)  (- — 0  <f>+(x)  dV  dVx 


p  x  m 


’  ii  p  x  m  p  n 


x  ()>*  (n)  4>->  (n)  dV  dV  dV 
n  n  p  T  n 


t ;  imn 


///  4>*  (p)  4>+(p)  ^  CT)C“~) 

i  P  r  m  pn 


x  (n)  (n)  dV  dV  dV 

r+v  ;  nv  ;  p  x  n 

n 


(II.  16) 


(II  Al) 


(II  .18) 


(II  .d9) 


the  integrals  P.  Q.  ■++ .  T,  and  Y,  +++  are  the 

6  l;im’  xt;im’  4- ;  imn  4  ;imn 

exchange  counterparts  of  the  above  integrals .  The  exchange 
integrals  to  the  first  two  are  given  by  interchanging  the 
electron  labels  of  the  second  and  fourth  orbitals.  For 
the  second  two  cases  the  integrals  are  obtained  by  inter¬ 
changing  the  electron  labels  of  the  fourth  and  sixth  orbitals. 

The  coefficients  appearing  in  Eqs .  (II. li)  to  (II.  14) 
are  constants  depending  on  the  occupancies  of  the  orbitals 


■ 


■ 

T; 

. 

. 
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involved  and  on  the  multiplicity  and  orbital  angular 
momentum  of  the  multiplet  under  consideration. 

i  2  i 

<$|H  I $>  may  now  be  written  as 


u2 

<H  > 


=  2  y  fXX-^  -  2  y  fXfy  +2  y  fXfy 

lii  L  lmimmi  L  l  m  l l  mm 

lm 


lm 
if  m 


+  y  fxfy(p^  +  q^)  -  y 

L  i  mv  im  xim  L 

-y-y  -y-y- 


fXfyfv 

l  m  n  lm  mmn 


ini 


lmn 
if  m 


+  y  fXfyfv(3H^N->->++  +  T-w  +  Y ->->->) 
^  i  m  n^  li  mmnn  lmn  lmir 

lmn 


y  f^fyfvf  c  (2N-^->-*N+-^) 
L  l  m  n  r>  nmn  nmpp' 


imnp 

m^n 


nmpp- 


y  fXfyfvfc{N^-^N 
L  i  m  n  p  limm 


->->  -  (  -  N  ■?—*"*•■*■[  2  N-* +  N-> ->->->  ]  )  } 
nnpp  ^3  imnp L  pmni  nmpijy 


imnp 


(11.20) 


Since  A  =  <H2>  -  E2,  E2  must  be  subtracted  from  (11.20) 
in  order  to  obtain  the  final  result. 


E2  =  4  I  fXfy  +  4  I  fXfyfv  H-J-+ 

. L  lm  n  mm  . L  lmn  n 


H-^N- 


lm 


mmnn 


lmn 


+  y  fXfyfvfc  N++++K++++ 
L _  imnp  limm  nnpp 

imnp 


(II.  21) 


After  subtraction  of  (11.21)  from  (11.20)  one  obtains, 
after  rearranging  terms, 


« 

I 
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A  =  2  y  f*X^  +  y  fXfy[P^->  +  0++  -  2H^+H->^]  +  ]  f*fyf 

L  ii  . L  l  m  lm  xim  lm  miJ  L  l  m 


-)• 

l 


lm 


lmn 


^  _p  y  4?  ^ 

n 


x  [ T +  Y ->-*->  ]  +  V  f  Af  y  f  f  ^ 

L  lmn  lmn  un  minnJ  l  m  n  p 

lmnp 


x  [=  N *— >■  (  2N'>->~>"}‘  +  N -*■■>■ ->•■*•)  -  2  N  ■►N  ■*•■*■■>■■*■  ] 

L3  lmnp  v  pmni  nmpi  nmn  nmppJ 


(11.22) 


If  now  the  term  in  the  first  bracket  is  denoted  as  A* 

* ,  lm* 

the  term  in  the  second  bracket  as  B*  and  the  one  in 

* ,  lmn 

the  third  as  C*  then  A  may  be  written  finally  as 

*,imnp*  J  } 


A  = 


2  7  f*X.  +  l  fXfyA*  ^  +  i 

L  ii  L  l  m  * , lm  L 

-V  -V  -V  7  -V-V- 


1 


lm 


fXfyfvB* 
l  m  n  x , lmn 


lmn 


+ 


l  fAfyfvfCC 

-i  m  n  t> 


lmnp 


l  m  n  p 


lmnp 


(11.23) 


If  one  now  defines  symmetrical  operators  with  respect  to 
subscripts  and  superscripts  [see  Appendix  II  of  the  paper 
by  Fraga  and  Birss  (29)]  the  form  of  Eq.  (2.6)  is  finally 
obtained.  It  must  be  noted  however  that  the  coefficients 
in  the  linear  combinations  of  integrals  in  (11.23)  and  Eq . 
(2.6)  differ.  This  difference  is  explained,  and  the 
relationship  between  the  two  sets  of  coefficients  is  to 
be  found,  in  the  paper  of  Fraga  and  Birss. 


' 


-10  2- 


Appendix  III 

In  this  Appendix  LCwdin’s  method  for  determining  lower 
bounds  will  be  presented  (  19  ) ;  his  notation  will  be  used 

throughout . 

It  is  necessary  in  order  to  understand  his  final  formula 
to  outline  its  derivation  here,  even  though  it  is  somewhat 
lengthy.  It  starts  with  the  partioning  of  the  secular 
equation  in  order  to  obtain  a  bracketing  theorem  for  the 
eigenvalues;  this  provides  a  means  of  giving  bounds  to  the 
eigenvalues .  Then  the  method  has  to  be  developed  in  terms 
of  an  operator  formalism  so  that  perturbation  theory  may 
be  applied  and  the  reaction  operator  introduced.  However 
the  reaction  operator  contains  an  infinite  expansion  and 
therefore  must  be  approximated  in  some  manner.  This  is 
accomplished  by  a  so-called  inner  projection  technique. 

With  this  accomplished  one  may  finally  write  his  formula 
in  its  final  form  after  a  number  of  algebraic  substitutions 
and  manipulations. 

Partitioning  the  Secular  Equation 

Consider  the  time  independent  SchrOdinger  equation 

H¥  =  ET  (in. 4) 

where  H  is  a  Hermitian  operator  and  y  is  the  exact  wave  func¬ 
tion.  Further,  suppose  E  is  the  lowest  eigenvalue  and  that 
it  is  discrete  and  non- degenerate .  Then  one  may  choose  a 
complete  orthonormal  set,  in  order  to  express  1  as 


■ 

' 


•  ’  f  '•  ■  V 
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?  8  I  CV4 


k  k 


(III. 2) 


Thus  equation  (III.l)  may  be  rewritten  in  matrix  notation  as 

(III  .3) 


H  c  =  E  c 


where 

Hki 

=  <$ 

klHIV- 

One  may  write  Eq 

.  (III. 3) 

then 

as 

M  c 

=  0 

(III. 4) 

where 

M 

=  H 

-  E  •  i. 

Now  suppose 

that 

the 

basis 

$  = 

{$^1  is 

divided  into 

two  sets, 

$  and  . 

a  b 

Then 

M  and  c  will 

become 

/Maa 

Mab\ 

-y 

/  C  \ 
a  \ 

ii 

2  « 

,  c  = 

\ 

(III  .5) 

'Mba 

Mbb' 

^  c  ; 

cb 

Hence  Eq.  (III. 4)  may  be  written  as 


M  c  +  M  i  c,  =  0 
.aa  a  „ab  b 


Mi  c  +  Mv  v  £,  =  0 . 

_ba  a  „bb  b 


(III. 6) 


If  one  solves  for  from  the  second  of  the  equations  (III. 6) 
then  ^  _  w-x  -  (III. 7) 


-*■  .  .-d  ..  ->■ 

c,  =  -Mi  ,  M,  c  . 
b  _bb  _ba  a 


This  equation  of  course  is  valid  only  if  the  matrix  M^  has 
an  inverse.  Combining  Eq .  (III. 7)  with  the  first  of  the 
equations  (III. 6),  one  obtains 

(Maa  -  Mab  Mbb  *Wc"a  " 


0. 


(III. 8) 
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Now  recalling  that  M  =  H  -  E  •  1  one  may  rewrite  Eq .  (III.  3) 
by  the  use  of  Eq.  (III. 8)  as 

l?aa  +  Hab<E  '  Ibb  '  Hbb^  Hba^a  =  E  Sa  Cm-9> 
The  subspaces  a  and  b  may  in  general  be  chosen  arbitrarily, 
but  for  the  present  application  the  choice  of  ,  i.  e., 

a  one  dimensional  subspace  for  a,  will  be  made.  With  this 
choice  Eq .  (III. 9)  becomes,  for  c^  f  0 

E  ”  H11  +  Hlb(E  •  !bb  -  Hbb5_1  Hbl-  c111-^ 

Since  this  gives  an  implicit  expression  for  the  energy,  E, 
LOwdin  defines  a  function  -  f(6)  given  explicitly  as 

*4  =  H11  +  Hlb«  •  hb  -  Hbb^d  Hbl  E  (hi .-Id) 

Equation  (III.lO)  is  reproduced  by  Equation  (III. 11)  when 
=  &  =  E.  Equation  (III. 11)  may  then  be  used  to  obtain 
a  "bracketing  theorem"  for  the  quantities  t  and  ti  . 

The  Bracketing  Theorem 

One  obtains  from  Eq .  (III. 11) 

d£j  - 

a.r  =  -Hlb<e  •  dbb  -  Hbb)  Hbi>  (i i 1 .12) 

and  from  Eq .  (III. 7)  one  has 

*  &  ■  dbb  -  Hbb^  Hbl  C1‘  (III  .13) 

Combining  Eqs .  (III. 12)  and  (III. 13)  one  obtains 

dti  =  _  ^b  ^b  .  (Ill  .14) 

d  t  c|  cd 


f.o  hr  ■  .- >■  :■  *af! :  m  .  adl  %n.i  t  -  ^s*rd 


■ 
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Now,  if  =  E  +  and  t  -  E  +  e  and  since  E  =  f(E)  one 
obtains  by  use  of  Lagrange's  mean  value  theorem 

E  +  =  f(E  +  e)  =  f (E)  +  e f '  (E  +  9e)  (111.45) 

hence 

=  e  f '  (E  +  6e)  =  e 

where  0  <_  9  4 .  Now  f'  (E  +  9e)  must  be  a  negative  quantity 

as  a  result  of  Eq .  (111.44),  hence  and  e  must  have  opposite 
signs  and  thus  and  t  must  "bracket"  at  least  one  true 
eigenvalue . 


-  jdf (E  j_gej  (III.l6) 

at  at 


The  Operator  Formalism 

In  order  to  proceed  it  is  necessary  to  write  Eq .  (III. 44) 
in  operator  form.  To  this  end  consider  the  projection 
operator  0  and  an  arbitrary  trial  function  $,  such  that 


0$  =  <j>  (III. 17) 

where  <f>  is  the  reference  function  which  spans  the  one  dimensional 
subspace  a,  i.  e.,  let  =  4> .  Let  the  operator  P  be  defined 
as 


p  =  4  -  0,  op  =  0 


where  d  is  the  identity  operator.  Then  (E 
be  considered  as  the  matrix  representation 
of  the  operator 

_4 

T  =  P[a  •  0  +  P(E  -  H)  P]  1 


Ibb 

in  the 


P 


(111.48) 
„  ,-d 

Hbb)  may 
basis  {t^}, 

(III  .19) 


where  a  j  0  is  an  arbitrary  number.  The  need  for  defining 


- 
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the  operator  T  in  such  a  manner  arises  from  the  fact  that 


(E  • 

i 

12 
12 
-v-i  ' 

-1  . 
l 

s  th 

e  inverse  of 

the  corner 

of 

a  matrix 

and 

it  must  be 

insured 

that 

thi 

Ls  inverse  properly 

exists . 

With 

the  defini 

t  ion 

of 

T  given 

above 

this  will 

be 

the  case , 

» 

Now 

by  defining 

an 

operator 

ft  = 

0 

+  THO 

(III, 

.20) 

such 

that 

m  = 

Eft 

ft2  = 

ft 

and 

Tr (ft)  =  i 

(III . 

.21) 

one 

obtains 

=  ft$ 

= 

<j>  +  TH  <j) 

(Ill . 

.22) 

since 

H¥ 

=  Hft  $ 

Eft$  = 

E'F  . 

(Ill  . 

.23) 

In  the  above  it  is  assumed  that  <$|0|$>  =  1.  Note  that 


(H  -  E)n  =  (0  +  P)  (H  -  E)n 

=  0  (H  -  E)ft  =  0  (H  -  E)  (0  +  THO) 


=  0 (H  +  HTH  -  E) 0  =  0 

which  gives  0E0  =  Q(H  +  HTH)0  (III. 24) 

If  this  equation  is  multiplied  on  the  right  by  $  and  on 
the  left  by  4*  and  integrated,  one  obtains 

E  =  <<J>  |H  +  HTH  |  <f>>  .  (Ill  .25) 


Hence  one  may  finally  rewrite  Eq.  (III. 11)  in  operator  form 


as 


tA  =  <<f>|H  +  H 


H|  <f>>  =  f  (£) 


t  -  H 


(III  .26) 


' 
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Introduction  of  the  Reaction  Operator 
Lbwdin  now  chooses  to  express  the  Hamiltonian  as 


v  H  =  H  +  V  (III. 27) 

where  Hq  is  the  unperturbed  Hamiltonian  and  V  is  the 

perturbation  and  the  reference  function  <J>  is  chosen  to 

be  d>  ,  then  with 

o 


one  obtains 


H  .  4>  =  E  <f) 

o  o  o  o 

0  =  I  <j>  >  <4> 

1  o  o 


(III  .28) 


TH<t>  =  THto  =  T(Hq  +  V)*o  =  T(Eq  +  V)  *0  =  TV*0 . (Ill . 29) 

Hence  Eq .  (III. 22)  becomes 

Tj  =  *0  +  TH  +  0  =  (4  +  TV)  4>o  (III. 30) 

and  Eq .  (III. 26)  becomes 

q  =  <*o|H  +  HTHUo>  =  Eo  +  <4>0|  V  +  VTV|*0>,  (III.31) 

and  t  =  V  +  VTV  is  called  the  "reaction  operator".  Thus 
Eq.  (III.31)  may  be  rewritten  as 


• 

A 

O 

-©- 

+-» 

O 

-0- 

V 

+ 

O 

W 

II 

<0 

(III .32) 

Since 

A 

T  =  P[a  •  0  +  P(C  -  H)  P  ]  XP  =  P[a  •  0  + 

PCfe 

A 

-  H  )P  -  PVP] 

(III .33) 

•  r _  j  _  • _ ^  _  nr  ^  ,  t./'A  tt  q 

A 

if  one  designates  P[a  •  0  +  P(&  -  H  )P]  P  as  T  , 

[a  •  0  +  P(6  -  H  )P]  as  A  and  [ P VP ]  as  B,  one  will  obtain 


the  following 

T  =  T  +  T  VT  =  T  (1  +  VT) 
00  o  J 


(III. 34) 


- 
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by  use  of  the  identity 


-1  -1  -1  -1 
(A  -  B)  1  =  A  1  +  A  (A  -  B) 


(III. 35) 


Multiplying  Eq.  (III. 34)  on  the  right  by  V,  one  obtains 

TV  =  T  (V  +  VTV)  =  T  t 

and  now  multiplying  on  the  left  by  V  and  adding  V  to  each 
side  of  the  equation 


or 


V  +  VTV  =  V  +  VT  t 

o 


t  =  V  +  VT  t 
o 


(III  .36) 


From  Eq .  (III. 36)  one  obtains 


t  -  VT  t  =  (1  -  VT  )t  =  V 
o  o 


and  dividing  both  sides  by  V  one  obtains 


d  =  (V"a  -  T0)t, 


-d 


provided  V  exists,  Hence  one  may  finally  write 

-1  -i 

t  =  (V  -  T  )  i. 

o 


(III  .37) 


However  T  involves  an  infinite  expansion  in  terms  of  the 
unperturbed  eigenfunctions  and  thus  t  cannot  be  determined 
explicitly.  Hence  one  must  try  to  approximate  t  in  some 
manner.  To  do  this  the  idea  of  an  "inner  projection" 
is  needed. 

The  inner  projection  of  t. 


Suppose  A  is  a  self-adjoint  pos it ive - definite  operator, 


then  one  may  write 


' 


% 


. 
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V 

A  =  I  akluk><ukl  (III. 38) 

ic 

where  a^,  and  are  the  eigenvalues  and  eigenfunctions  of 
the  operator  A,  and  a^  >  0 .  One  may  now  define  the  positive 
square  root  of  A  as 


Ad/2  =  l  ak/2|Uk><Ukl •  (III. 39) 

k 

Now  using  the  projection  operator,  0,  one  may  define  an 
"inner  projection",  A',  as 


A'  =  aM2OA^2 


(III .40) 


where 


0  =  fA  ift  =  I  l£k>Aki<f; 


(III .41) 


for  a  set  of  linearly  independent  elements  {f.},  and  A^  = 
<^k^l>kl*  £nner  Projection  is  such  that  it  satisfies 


the  inequality 


0  <  A’  <  A. 


(Ill  .42) 


Thus  if  t  is  a  positive  definite  operator 

t»  =  t^/2fA'1ftt^-/2  <  t  (III. 43) 

and  hence 

E  +<(j>  1 1 1  4  >  >  E  +  <  <}>  1 1 T  |  >  .  (Ill, 44) 

o  o 1  1  o  o  o 1  1  o 

Thus  t  ’  will  provide  a  lower  bound  to  t  and  6^  =  E  +  <<t>  |t'|<t>  > 
will  provide  a  lower  bound  to  the  true  eigenvalues.  If 
one  substitutes  h  -  t1/2i  into  Eq .  (III. 43)  one  obtains 

tf  =  hA  hT,  A  =  hTt  1h .  (III. 45) 


Lowdin  then  chooses  the  substitution  h  =  (fe  -  H  )t,  which 

o  j 

gives 


. 
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?  -A 

<  <J)  t '  I  4>  >  =  (t  -  E  )  <  4>  |  i  >  A  I  <f)  >, 

o  1  o  o  oIJ  Jlo 


(III. 46) 


where 


j  -  a4.  j  ^ ,  .  ..,  jR)  is  an  arbitrary  function  space 
-1 


and  where  A  is  the  inverse  of  the  following  nth  order 


matrix : 


A  =  <t  |  (6  -  Ho)Vd(c‘  -  Hq)  |  j  >  -  <j  1 4  -  Ho  I  j> 


+  (6  -  E  )  <T  I  <j)  >  <  t  |1>. 
v  o  1  o  oIJ 


(III  .47) 


V 


-1 


Application  to  the  Ground  State 

For  the  case  of  the  application  to  the  ground  state 
=  r^2  and  Eq  •  (III. 47)  may  be  rewritten  as 


A  =  «2<j|rl2lT>  -  fc<j|Hor12|T>  -  £<j|rl2Ho|T>  + 

-l  <Tl?>+  <1  |H0  |  J  >  +  [t  -  E0)<3U0><*0|j>  (III. 48) 
where  the  are  of  the  Hylleraas  type,  i.  e., 


r  •  i  r  •»  f\T  1  .  ^  1*1  “  S  f  2  i 

{Jk}  =  {xlmn}  =  {Nlmns  1  u  e  }' 


(111.49) 


The  evaluation  of  the  integrals  in  Eq .  (III. 48)  using  the 
functions  of  Eq.  (III. 49)  is  discussed  in  Appendix  IV. 

In  order  to  carry  out  a  calculation  of  this  type, 

one  must  first  compute  the  matrix  A,  then  find  its  inverse, 

-i 

A  ,  one  must  then  carry  out  the  multiplication  of  Eq .  (III. 46) 
and  finally  add  the  latter  quantity  to  the  unperturbed 
energy  in  order  to  obtain  the  lower  bound 


' 

' 


I 


-  Ill  - 


Appendix  IV 

In  this  Appendix  the  explicit  forms  of  the  integrals 
needed  in  the  application  of  the  method  proposed  in  this 
thesis  and  in  the  application  of  LOwdin’s  method  to  the 
groundstate  of  He  will  be  presented.  These  integrals  are 
evaluated  by  a  modification  of  the  analysis  of  Wilets  and 
Cherry  (36) . 

The  functions  used  in  these  calculations  are  of  the 


Hylleraas  type  and  are  given  by 

XT  l.m  n  -s/2 

y =  N,  s  t  u  e 
Almn  Imn 


(IV. 1) 


where  s  =  k  (r^  +  ,  t  =  k (r^  - . r^)  and  u  =  kr^ »  m  is 

always  an  even  integer.  It  is  convenient  for  what  follows 
to  define  the  following  functions 


A(p,  q,  r)  = 


2  7T  I*00 

~T6  J  0 


s  r>-  'u  tif.2  -  f2 


/"  e‘sspds  /Q  urdu  /“  tq( 


t  )  dt 


=  B(p  +  2,  q,  r)  -  B  (p ,  q  +  2 ,  r) 


B(p,  q,  r)  =  /“  e"SsPds  urdu  tqdt 

k 


(IV. 2) 


■  7T 


(q  +  4)k 

„  2 


In  e'SsPds  /*  uq  +  r  +  < 


6  J  0 


du 


■  TT 


(q  +  i) (q  +  r  +  2)k 

=  2u2(p  +  q  +  r  +  2)! 
k^  (q  +  1)  (q  +  r  +  2) 


in  e  s 


~  s  p  +  q  +  r  +  2 


ds 


(IV. 3) 


1 


* 
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C(p,  q,  r)  =  /“  e " SsPds  JSQ  u*du  /“  - -j  dt 


s  -  t 


However , 


/ 


u  t  1 
0  s2  -  t2 


q/2  -1  2i  0  -  a  , . 

dt  -  -  l  -^-[t^2l-V  +  /“ 

i  =  0  q -  2 i -d  U  s  -  t 


q/2-d  2i  q-2i-l  v 

n  "  s  un  ,  q  -  1.,  i  s  +  u 

-  +  s n 


I 

i  =  0  q-2i-l 


In 


s  -  u 


,q  q/2-d  2i  0.  j 

r s  r,  ru  tn  j.  n  v  s  f s  q+r-2i-l 

/n  u  du  in  — - T  dt  =  -  l  - —  /  u4 

1  =  0  q - 2i  -1 


0 


0  s2  -  t2 


du 


q  -1 

,  rs  r  -i  |  s  +  u  |  j 

+  -  J0  u  In  I  -  7  —  \  du 


s  -  u 


=  -s 


q  +  r 


q/2  -1 

l 


d 


sq  +  r  Xt  id  +  z 


i  =  0  (q-2i-d)  (q  +  r-2i)  2 


+  5 -  /£  zUn 


4  - 


dz 


Hence , 

C(p,  q,  r) 


2 

"5"  J  0 


q/2-d 


2„  sp*q+re-sds[.  -  j 


i=0  (q-2i-l) (q+r-2i) 


1  1 1  1 1  +  z  |  j  t 

+  2  <o  z  ln'r^T ldz]  • 


However,  because  the  Hylleraas  expansion  is  used,  q  is 
always  an  even  integer  and  greater  than  zero,  hence 


q/rd  4 

i  =  0  (q-2i-l)  (q  +  r-2i) 

d 


d 


q/2-l 

I  — 

j  =0  (2 j  +1] (2 j  +  r+2) 


k=2  (k-1) (k+r) 

k , even 


ul 
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The  logarithmic  integral  is  treated  in  the  following 


manner 


DO)  =  \  Iq  z  r  In  \\  ~j— ~  1  d  z  ■  =  |  z  r  In  (1  +  z  )  dz  -  |  Jq  zrln(d-z)dz 


=  ^  /2  (x-l)rlnx  dx-i-^J —  A  (x-l)rln  x  dx 

Z  i  o  0 


2 

=  1  I  (£) ('l)k [ (-1) r  /?  Xkln  x  dx-/^xkln  x  dx] 
Z  k=0  10 


1 


k+1 


2  l  (^)(-d)K{(-l)r[^— m  x-  xk-YT]d 

z  k=0  K  k+1  (k+l)z  1 

k+1  k+1  a 

r  X  ,  v  X  nln  , 

-  [ — —In  x - — T]  ]  } 


k+1 


(k+1)2  0 
,  k+1 


1  r  ,rv  ,  ^kr  ,  jsrr2  ,  c 


k  +  1 


2  l  (1)(-1)"{(-1)X[ 

z  k=0  K 


k+1 


-In  z  - 


1 


(k+1)2  (k+1)2 


] 


1 


(k+1) 


\  i.  ♦  i^o[2k+hn  z 


k=0 


(k+1) 


k+1 


2k+l  4 

i -  +  -L_] } 

k+1  k+1 


for  r  f  -1,  for  the  case  r  =  1,  then 


D(r)  =  D ( -1)  = 


n 


8 


Thus,  finally  C (p ,  q,  r)  may  be  written  as 


c(p,  q,  r)  =  2-7T-  Cp+q+rII.{-  ^  - A -  +  D(r)}  (IV. 4) 


k=  2  (k-1) (k+r) 

k  even 
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All  of  the  necessary  integrals  may  now  be  evaluated  in 
terms  of  these  functions . 


.  i  -  AT  .  jw  2k  -v  l+m+n+6  fT,T 

<<t,o'xlmn>  =  (~J  Ni™«AC1.  m>  n+d)( - )  (IV.  5) 


k  +  4 


where 


N,  =  [ 
lmn  L 


d 


A(l,  m,  n+d) 


d 

2 


<f>  =  ^-e  2s^  for  the  case  of  the  ground  state. 


<xlmn  'T/l2'xl,m,n,> 


N,  N, ,  ,  . 

lmn  1 1 m 1 n ' 

~~  kP 


A(L,  M,  N  +  p) 


where  L=l+l,,M=m+m,  and  N  =  n  +  n*  +  d;  two  cases 
exist,  p  =  0  and  p  =  d. 

For  the  He  atom  the  Hamiltonian  may  be  written  as 


.2  . 2  -  2  n  ,  2  .2.  »  2  2. 

H  =  -k2[-^-  +  +  2s(u^  -  tj  9  +  2t(s  -  u  )  3 


3s2  9t 2  3u2  u(s2  -  t2)  3  s  9  u  u(s2  -  t2)  3 1 3u 


4s 


4t  L  +  1  L  +  _  4Zs 


?  ?  7  7  72 

-  t  3 s  s  -  t  3t  u  3u  k(s^  -  t  ) 


]  + 


u 


If  the  term  in  brackets  is  denoted  as  H',  then  -k  Hf  =  H 
where  H^  is  the  unperturbed  Hamiltonian.  Making  use  of 
the  following  equations 


0 
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_a_  =  ri  _  4. 

3s  xlmn  As  2^xlmn 


3 


:  2  lmn 
3s 


|i(i  -  i)  .  i  +  iu 

L  2  s  4 J  xlmn 

s 


3  _  m 

3T  xlmn  t  xlmn 


m(m  -  4) 


3 1 


2  Almn 


■lmn 


3  n 

3u  xlmn  u  xlmn 


2  xlmn 


n  (n  +  4) 


3u 


u 


lmn 


1  4,n 


o  _  /- j.  x-\jli  a  _  mn 

3s  3u  xlmn  ^s"  2Ju  xlmn  3t3u  xlmn  7u  xlmn 


one  may  write 


Hlv  >  =  ri<2(I  1C1  -  4)  i  m(m  -  4)  n (n  +  4) 

|xlmn  LK  *s  2  4  ^2  "  2 

s  t  u 


[Cj  +  1  -  *f)s  -  (1-m)  (n+2)  +  nl  H 

u 


2  ,  2 
s  - 1 


.  2  2  , 
n  st  s  -j  i  ,  k -j 

-  mn  -j]  }  +  -]  x 


2  2 

u 


u 


uJ  Almn 


k  xlmn{1[''1>0*0]  '  |[°>0>°]  ‘  W  -  'D  [-2,0,0] 


-  m(m  -  1)  [0,-2 ,0]  -  n(n  +  d)  [0,0, -2]  +  1[0,0,-1] 


+  -riKn  +  2  -  Hi  [4,0,0]  -  2(1  -  m)  (n  +  2)  [0,0,0] 

S  -t  k 


+  2nl[0,2,-2]  -  n  [4 , 2 , -  2 ]  -  2mn [ 2 , 0 , - 2]  ]  } 


where  [a,  3,  y]  =  sat^uy. 


' 


V 
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One  may  now  write 

<  x  1  mn  ^  rl  2 ^  x  1 '  m ’  n  ’  > 

x  A(L  +  aF,  M  +  3F, 

x  B (L  +  a?,  M  +  8?, 

<xlmn  ^orl2  ^  X1  'm’n  ’  > 
x  A(L  +  a?,  M  +  8?, 

x  B(L  +  a*? ,  M+  3*? , 

N 

<xlmn  ^  ^  X1 ’ m'n ’ > 
x  A(L  +  a?,  M  +  8?, 

x  B(L  +  a?,  M  +  8?, 


' m ' n ' ^lmn 


N  +  y?  +  1)  + 


N  +  +  1) 


i  Ni 

lmn  1 
k 


'm’n 

5 


{ 


N  +  +  1)  + 

N  +  y^  +  1) ) 


i  Ni 

lmn  1 


k 


'm'n' ^ 
4 


N  +  Y^)  } 


5 

T  F-Cl'm* 
i=l  1 

5 

l  G. (1 ' ,mf , 
i  =  l  1 


F±  Cl ' ,m» ,n ’ ) 
Gi(l,m'n') 


<Xt  |Hr,9H|x1t  r  i  >  ~ 
Almn 1  o  12  o'l’m'n’ 

5 

[  l  F.(l!m;n'+1)A(L 

3=1  3 


Nn  N-,  ,  f  ,  5 

lmn  1 ’ m ’ n ’ r  r 


{  y  f- ci, 

i=i  1 


»  t 

m. 


+  aF  +  ctF,  M  +  gF  +  3F, 

i  j  ’  i  y 


n') 

(IV. 7) 

n’  +  1) 

’  +  1) 

(IV.  8) 

(IV. 9) 

nVl) 

N  +  yF  +  yF  +1) 
i  1 


X 


- 


r  ' 

' 
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5 

J, 

5 


+  Gj (1 ’ mj n ' + 1 ) B (L  +  cu  +  aV,  M  +  B?  +  3^,  N  +  y^  +  yV  +  1)] 

5 


+  l  G-  Cl’.m'.nVl)  [  l  F  .  (1 JmJn 1 +1 ) B  (L  +  a  +  a 

i=4  j=l  J 


G  .  F 

j 


M  +  3?  +  3 j  ,  N  +  y?  +  yF  +  1)  +  ^  Gj  (1,'mJnJ  +  l) 


C (L  +  a?  +  a?,  M  +  3^  +  3?,  N  +  y^  +  Yj  +  1)])  (IV. It)) 


<Hxlmn|Hxl 


'  m '  n  ' 


>  = 


N,  N-.  ,  ,  ,  6  6 

-lmn  L  Fi(l,m,n)[  l  F  (1 '  ,m’  ,n ' ) 

1  •  A 


A(L  + 

a?  + 
l 

F 

V 

M  + 

■Z- 

SF 

V 

N  + 

B(L  + 

aF  + 

l 

G 

V 

M  + 

bg 

V 

N  + 

rl' 

F-. 


1=1 
5 

L 

j=i 

5 


[  F.  (1*  ,m»  ,n')B(L  +  ai  +  a  j  >  M  +  3?  +  3?,  N  +  y?  +  yF) 

+  ^  Gj (1 ' ,m' ,n ’ ) C(L  +  a?  +  a?,  M  +  3?  +  3j ,  N  +  y^  +  y? ) ] } 

(IV. 11) 


where  the  above  quantities  a^}  3^,  y^,  F^  and  G^  are  given 
in  the  table  below.  It  should  be  noted  that  the  value  of 
the  integral  <X£mn | Hq | X£ * m » n ' >  is  found  by  using  Eq . 

(IV. 9)  but  using  the  limit  of  the  first  summation  as  5_ 
rather  than  6. 


, 


■ 


Table  4.2 


i 

F 
a  . 
l 

F 

Yi 

F^ ( lmn) 

G 
a  . 
l 

rG 

Bi 

G 

Yi 

G^ ( lmn) 

1 

0 

0 

0 

~4 

0 

0 

0 

-  2 (1 -m) (n+2) 

2 

-1 

0 

0 

1 

1 

0 

0 

(n+2  -  4Z/k) 

3 

-2 

0 

0 

-id-1) 

0 

2 

-2 

2  In 

4 

0 

-2 

0 

-m(m-l) 

1 

2 

-2 

-n 

5 

0 

0 

-2 

-n (n+1) 

2 

0 

-2 

-  2mn 

6 

0 

0 

-1 

l/k 

With 

the 

expressions 

given 

above 

all 

the  integrals 

needed  for 

the 

two 

methods 

may  be  computed 

• 

